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CALCULAT IV E TECHNIQUES. FOR TRANSONIC FLOWS 
ABOUT CERTAIN CLASSES OF WING-BODY COMBINATIONS — PHASE II 

By Stephen S. Stahara and John R. Spreiter* 

Nielsen Engineering & Research, Inc. 

SUMMARY , 

Theoretical analysis and the development of associated computer 
programs were carried out for the purpose of developing computational 
techniques for predicting properties of transonic flows about certain 
classes of wing-body combinations. The procedures used are based on 
the transonic equivalence rule and employ either an arbitrarily-specified 
solution or the local linearization method for determining the nonlifting 
transonic flow about the equivalent body'. Theoretical results obtained 
by using the local linearization method are presented for surface and 
flow-field pressure distributions for certain members of the general 
classes of configuration's studied, for both nonlifting and lifting 
situations, at M^ =1. 

The computational 'programs developed under this report are documented 
and presented in a general user's manual included as part of the 
report. 

INTRODUCTION 


Stimulated by the need for accurate prediction of transonic flows 
about realistic aircraft configurations, recent research is producing 
significant advances in the ability to predict theoretically both two 
and three-dimensional transonic flows about a wide variety of aerodynamic 
shapes. While current emphasis seems to be placed on the development 
of numerical techniques (refs. 1, 2, 3, 4, 5) , it has become clear that, 
although significant accuracy limitations need not exist for advanced 
computer programs, these techniques do have cost limitations with regard 
to both accuracy and the use of alternate methods. Consequently, in 

*Professor, Departments of Applied Mechanics and Aeronautics and Astro- 
nautics, Stanford University, Stanford, California. (Consultant at 
Nielsen Engineering & Research, Inc.) 



order to enhance these computational efforts, the parallel development 
of proven analytical and analytic/numeric methods to provide accurate 
first approximations, for example, in the systematic study of a large 
number of configurations, is clearly warranted. 

Previous investigations by Nielsen Engineering & Research, Inc. 

(NEAR) in reference 6, where the local linearization method and the 
transonic equivalence rule were applied to predict surface and flow-field 
properties of several general classes of axisymmetric and nonaxisymmetric 
bodies for both lifting and nonlifting situations, and in reference 7, 
where those results were extended to include several classes of wing- 
body combinations, have demonstrated the effectiveness of such a com- 
bined approach . 

While the ultimate goal of the present investigation is to develop 
computational techniques for the prediction of the flow field, pressure 
distribution, and aerodynamic characteristics of three-dimensional, lift- 
ing, wing-body combinations, the purposes of this study are to extend 
the results of reference 7 to include a more general class of wing 
planform shapes, specifically, wings haying (1) sweptback trailing 
edges, and (2) finite tip chords. In addition, the computer programs 
developed in reference 7 were to be further enhanced with regard to 
minimization of computational time and applicability to wider classes 
of equivalent body shapes and equivalent body transonic solutions. 
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lateral distance from wing leading edge 
wing ordinates, eqs. (8), (9) 

angle of attack 

wing leading-edge sweep angle 

wing trailing-edge sweep angle 
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polar angle in crossflow plane 

ratio of major to minor axes of elliptic cross 
section, a/b 
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thickness-to-chord ratio of wing profile, eq. (10) 
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perturbation velocity potential associated with two- 
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and growth of cross section in crossflow plane 

perturbation velocity potential associated with two- 
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contraction of equivalent cross section in cross- 
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perturbation velocity potential associated with two- 
dimensional incompressible solution for expansion 
or contraction of cross section in crossflow plane 

perturbation velocity potential associated with twprr 
dimensional incompressible solution for translation 
of cross section in crossflow plane 



ANALYSIS 


General Considerations 


Because the current work is an extension of that of reference 7, 
the basic theory and equations used are discussed in depth in that 
reference and their derivation will not be repeated here. For convenience, 
however, those points relevant to the present work will be outlined. 

The coordinate system used for all of the three-dimensional flows 
considered herein is a body-fixed Cartesian system centered at the 
body nose with the x axis directed rearward and aligned with the 
longitudinal axis of the body, the y axis directed to the right, 
facing forward, and the z axis directed vertically upward, as shown 
in figure 1. For lifting situations, the free-stream direction is 
taken to be inclined at any arbitrary small angle a to the x axis 
and confined to the' x-z plane so that there is no sideslip. The 
governing partial differential equation for the perturbation potential 
0 is given by t 

M m 2 (y + 1) 

(1 - M 2 )0 +0 +0 = — 0 0 ' ( 1 ) 

oo ^ XX ^yy ^ ZZ X XX 


where M is the free-stream Mach number, y 
heats, and the free-stream velocity. The 

Cp in the above reference frame is given by 


the ratio 
pressure 


of specific 
coefficient 
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(^y 2 + 0 Z 2 ) 
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The transonic equivalence rule enables the perturbation potential 0 to 
be expressed in the form 


0 - 02, a + 02, t “ 02, B + .0 


(.3.) 


where each of the individual components has the meaning indicated in 
figure 1. Since 0 2 , a , 02 t> an< 3 02, B satisfy the two dimensional, 
incompressible Laplace equation 
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yy ZZ 

(where the subscript l depends upon the particular potential in 
question), they are independent of Mach number. Hence, the only portion 
of the solution dependent upon is and this term represents 

the solution to the full transonic equation (1) for the nonlifting flow 
about the equivalent body of revolution. Because the equivalence rule 
places no essential restrictions on the methods of calculating 0 0 , its 
solution may be determined in a variety of ways. For example, it can 
be given by an exact numerical solution, by experimental data, by an 
approximate analytic solution, or by a combined analytic/numeric solution 
such as the local linearization method. One of the tasks of the 
present work is to extend 'the applicability of the computer programs 
developed in reference 7 to include general, arbitrarily-specified 
solutions for <p , and the method of doing this is detailed in the 

user's manual. If the local linearization method is used to determine 

\ 

, or more conveniently, u D = (0„) , then one of the following set of 

B B B x 

three first-order nonlinear differential equations must be integrated 

according to whether ~ 1 , < M cr or M cr u < M^, where M 

M are the lower and upper critical Mach numbers, respectively. 
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and for purely supersonic flows (M < M ) 

^ 1 c cr,u oo 


d 

dx 



S'", (x) 
_£b 

4tt 


ln 


(M 2 
X 00 


1 + k V 


d 

dx 


r g ; b (x) 

47 r 


ln 


S ab (x| 

47TX 2 


1 

2tt 


A 

/ 


eb 


(x) - 




x - ? 


d? 


(7) 


where C in equation (5) is Euler's constant ~0.577, k in equations (6) 
and (7) is equal to M^ 2 (7 + 1) /u m , S ^(x) represents the area distri- 
bution of the equivalent body, and primes indicate differentiation with 
respect to the appropriate variable. These differential equations have 
been programmed for computation in reference 6 where details regarding 
starting conditions, numerical techniques, accuracy, limitations, etc. 
are provided. 
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Wing and Body Geometry 

The classes of wing-body configurations examined in reference 7 
and in this study are composed of finite thickness wing and either 
circular or elliptic cross-sectional bodies in which the bodies are 
area-rule indented along the wing-body junction in such a manner that 
the total cross-sectional area distribution (body plus wing) remains 
identical to that of a smooth body having a certain specified profile. 

The general class of wings considered have symmetric planforms consisting 
of straight leading and trailing edges, swept at arbitrary angles (3 

x 6 

and 0 respectively, to the y axis. In reference 7, the planform 
shapes were restricted to wings having either straight or sweptforward 
trailing edges and zero taper ratio. This work extends that class to 
wings with sweptback trailing edges and taper ratio between zero and 
one, as shown in figures 2 and 3. The wing profiles are represented 
by expressions of the form 


or 



( 8 ) 



t m (m/m - l) 
w 

2 (m - 1) 



x 

c 

w 



(9) 


where c y is the local chord, x the distance from the leading edge, 
m is a constant ^ 2, and t w is the wing thickness-to-chord ratio. 

In addition, the wings are assumed to maintain a constant thickness-to- 
chord ratio across the span, with the consequence that the wing profiles 
at all spanwise locations are geometrically similar. Thus, 


w 
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(Z (x,y) ) 
w ’ J max 

c (y) 
w Jr 


(Z w (x > Q)) max 

C R. 


where C Rt is the wing root chord. 


( 10 ) 
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Two categories of body shape are considered. Figure 2(a) and (b) , 
illustrates two members of the first category which have indented bodies 
that are circular in cross section, while figure 3(a) and (b) illustrates 
two members of the second category which have indented bodies that are 
elliptic in cross section and that maintain a constant ratio A(=a/b) of 
semimajor to semiminor axis along the entire body length. In reference 7, 
the profiles of the equivalent bodies of revolution of the wing-circular 
body combinations are described by the expressions 



n/n (n-i ) 

T . n 

eb 

2 (n - 1) 



( 11 ) 


where the exponent n is given in terms of the location of maximum 
radius by 
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where 


(ft 


1 / (n-i ) 


max 


(14) 


while the equivalent bodies of the wing-elliptic body combinations are 
parabolic-arc bodies, i.e. equations (11) or (13) with n = 2. This 
work extends the class of equivalent body profiles for both the circular 
and elliptic body shapes to include arbitrarily specified functions 
subject to certain closure and continuity restrictions on the derivatives 
that are discussed in the appropriate section of the included user's 
manual . 
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• Straight or Sweptforward Trailing' 

. Edge Planforms 

Circular bodies .- For finite thickness wing-circular body combina- 
tions having wings with finite tip chords and trailing edges that are 
either straight or sweptforward, the complex potentials, W 2>a , W P , and 
W 2 , b can be readily determined from the work of reference 7. 
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where 


is the complex variable in the crossflow plane 


o = y + iz 


(18) 


R^ is the indented body radius, s - depending upon axial location - 
represents the local wing semispan of either the leading (s = s^) or 
trailing (s = s fc ) edge, and represents the wing ordinates. The 
velocity components associated with these potentials can be found by 
substituting those expressions into the general formulas: 
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where the subscript i depends upon the particular potential in question 
and R.P. signifies the real part of a complex quantity. 

These operations have been carried out and the resulting expressions, 
which are quite lengthy, are given in reference 7. It should be noted 
that, in the evaluation of the velocity components associated with the 
thickness problem (W 2 ^-) , different sets of expressions are necessary 
depending on whether the point of interest is (1) at a general location, 
(2) on the wing surface, or (3) at the wing-body junction. These 
distinctions are required in order to account properly for the Cauchy 
singularities which appear in several of the integrals associated with 
the thickness velocity components. No such distinctions are required 
for the lifting (W 2jCl ) or equivalent thickness (W 2>B ) problems. 

Because of the symmetry of the class of wing-body combinations 
considered, nonlifting flows will produce no lateral forces or moments. 

The only force will be the longitudinal drag force which can be 
determined through the general formula, 


D 


a=o 



i 


$2 , B 



( 21 ) 


where D ^ represents the drag of the equivalent body while the other 
two terms involve the line integral along their respective contours 
(Cj_ is the contour defined by the cross section of the wing-body 
combination while C D is the contour about the equivalent area circular 

n , 

cylinder) of the product of the appropriate velocity potential and the 
normal velocity associated with it. We note that the drag indicated 
by equation (21) refers to the inviscid drag of the configuration minus 
the base pressure drag. As pointed out in reference 8, there exist 
many shapes of aerodynamic interest for which the two integrals involved 
cancel. In particular, we note that if the equivalent body and the 
original body have the same shape and surface slope at the base, as is 
the case for configurations studied here, then since both integrals 
are carried out over the same contour along which 0 2 t = 0 2 , B anc * 
d0 2 , t/5n = S0 2j B/dn, ;the integrals cancel and D a=0 = D e b* we 

define a drag coefficient C D based upon the maximum cross-sectional 
area of the equivalent body, S , then 
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where is the axial location of the body base and Cp e ^ is the 

pressure coefficient on the surface of the nonlifting equivalent body 
and is equal to 



(23) 


For the lifting situation, an exact analysis of the aerodynamic 
forces and moments, even within the framework of small disturbance 
theory, is quite formidable. The general formulas for determining the 
coefficients of lift, pitching moment, and drag are given by 
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where now the contour C, while still taken at the base of the body, 
must now account for the vortex wake which springs from the wing trailing 
edge and, as before, the drag given by equation (26) represents the 
inviscid drag minus the base pressure drag. Because the vortex lines 
near the body surface must follow the streamlines of the flow around 
the body, the vortex wake will not proceed parallel to the x axis, 
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in general, as it does in many simpler cases considered in slender body 
theory; but will move away from or toward the body axis to follow the 
lateral expansion or contraction of the flow field near the body as shown 
below. 

r 



The resulting flow at the body base is influenced by the wake and, 
consequently, is no longer independent of the flow at cross sections 
preceding it. The solution of problems of this type is discussed 
briefly in reference 9. In general, they are quite difficult to solve 
and since they are by no means unique to transonic slender body flows, 
their exact solution is clearly beyond the scope of the present investi- 
gation. Because the analysis presented here, however, remains valid up 
to the axial location of the wing tip trailing edge x = ^ sm2 (i.e. as 
long as the edge of the wing remains a leading edge) an estimate can be 
made of these quantities by making the assumption that beyond that point 
the vortex sheet remains parallel to the x axis and does not vary 
with x. With this premise in mind, we can proceed to evaluate 
equations (24), (25), and (26). Carrying through the indicated operations 

(see ref. 7 for details) , we arrive at the result that the coefficients 
of lift, drag, and pitching moment are given by 
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where the drag coefficient at zero lift c D a=0 is 9 iven 'by equation (22) 
and x r^ei < as ax i a l location of the point where the wing leading 

edge pierces the body surface. 

Elliptic bodies .- The basic analysis of wing-body combinations 
composed of wings having finite tip chords with trailing edges that are 
either straight or sweptforward and bodies having indented elliptic 
cross sections such that the total cross-sectional area distribution equals 
the area of the original smooth body with elliptic cross section proceeds 
in a manner analogous to that used for the circular body shapes. 

Apparently, the most direct approach consists of reducing the elliptic 
cross section to a circular one by use of the appropriate Joukowski 
transformation and then applying the methods used for the circular 
shapes. The transformation required to take the ellipse into the circle 
shown below 



is given by 
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Rx = 


(a + b) 


and the semispan ,s into the shortened semispan 


s + V s 2 - c 2 
2 


The potentials W 2jC[ , W 2 j t, and W ?jB are then given by 


2 > a ■ f x 

11 — = “ ia < °i + TT" 


(' ' ¥)']' 


_ii = i r 

'oo " *' 1 


- Sl + 4 T 72 ! 


(c 1 2 - ^ 1 2 ) ( o i ^ 2 ) 


+ 2 V x ,a> ||) i" o. 


w c,B S Ab (x) 


In o 


The velocity components associated with these potentials are found 
through the operations ■ 
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These have been carried out and are given in reference 7 where it must 
be remembered in those formulas that depending upon axial location s-l 
represents either the leading (Sj = s 2 ) or trailing (s x = s lt ) edge 
in the transformed plane. For the case of nonlifting flows about these 
classes of symmetric configurations, no lateral forces or moments exist 
so that the only force present is the longitudinal drag force. This can 
be determined through the use of equation (21) where now the contri- 
butions of the two line integrals do not cancel since the contour over 
which the product 0., i- c)0 2 i-/Sn is evaluated is the elliptic cross 
section at the base of the body whereas the contour for evaluating 
0 2 ,b is t,ie circular cross section of the equivalent body. 

Carrying out the indicated operations, we find that the drag coefficient 
of this general class of nonlifting elliptic wing-body combinations is 



(40) 


where Cp ^ i- s the drag coefficient of the nonlifting equivalent body 
and is given by equation (22), and K(FJ is the complete elliptic 
integral of the first kind. 

For lifting flows at small angles of attack about these configura- 
tions, if we apply the same assumptions regarding the trailing vortex 
wake as were, made for the circular body case, then the evaluation of 
equations (24), (25), and (26) provides the following results for the 
lift, pitching moment, and drag coefficients. 
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where c D a=0 i- s t ^ ie drag coefficient at zero lift and is given by- 
equation (40) . 
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Sweptback Trailing Edge Planforms 


Circular bodies .- For finite thickness wing-indented circular 
body combinations having sweptback trailing edges, as illustrated 
below, 



new potential solutions are required to account for the gap between 
wing and body cross sections which appear in the crossflow plane between 
the axial locations x = X rte x and x = x sm 2 * For lifting flows, the 
analysis is complicated even with the previously made assumption that 
the vortex sheet eminating from the trailing edge remains parallel to 
the x axis beyond the point x = X smg . This is due to the presence 
of the vortex sheet in the gap between the wing and body from X rt6i to 
X S m 2 - Consequently, i the flow at any axial station beyond X r ^ ei is 
influenced by the wake ahead of it and is no longer independent of the 
flow at preceding cross sections. Thus, the simplified analysis which 
was valid for lifting flows about wings with straight or sweptforward 
trailing edges - i.e. for cases where the edge of the wing always 
remained a leading edge - does not apply here. A new potential solution 
for 1 ^2 } d ' ' is- ; required and this is beyond the work scope of the present 
investigation. •• . 


r ( 
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One of the primary goals of this study, however, is to determine 
the potential w 2 t associated with the thickness problem for configura- 
tions of this type. For x < X rt6i , equation (16) is valid. For 
x rtei < x < X sm ^ , a new potential must be developed to account for the... 
gap between the wing and body. This is accomplished by using an extension 
of the method developed by Stocker in reference 10. That method is 
based upon the method of singularities and models the wing thickness by 
placing a continuous distribution of two-dimensional incompressible 
sources (or sinks) along the wing chordal plane together with their 
appropriate images within the body. The body itself, is represented by 
a source (or sink) at the origin. Although originally developed for a 
wing attached to a body, this method can accommodate such a wing body 
as shown below 



by distributing sources (or sinks) only along the wing and appropriately 
imaging them within the body. This method thus provides the following 
expression for the complex potential w 2 ,t( x >Y> z ) : 



Since for sweptback trailing edges with X rtei < x < x sm 2 > 
equivalent body area distribution and actual body area distribution are 
related through the expression 
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The velocity components associated with this potential are deter- 
mined through equations (19) and (20). As in the case of straight or 
sweptforward trailing edges, proper account must be taken of the 
familiar Cauchy singularity which appears in several of the integrals 
associated with these velocity components on the wing surface. If, 
because of the symmetry of these configurations, we restrict attention 
to the first quadrant of the crossflow plane (y ^ 0, z ^ 0) , the 
following results are obtained. For points at any general location but 
not on the wing surface, 


22 





23 




rf dZ w (x,P). 


z 2 + ( y - f, ) 2 z 2 + ( y + ■'; 


rT 2 

z 2 + (y - — : — ) 2 


: . ... V. . ■ 

. y + — 

+ + 

Z 2 + (y + — S_) 2 


z 2 + y 2 


, s 


%dZ w (x,^) S eb (x) 

dr + 

dx J 2 


Vt _ 1 f f* d2 w (x »^ ' 

7r [ Z J d * 

'» Sx. 


1 1 

+ 

z 2 + (y - z) 2 • z 2 + (y + e, ) 


Z 2 + 


R, 2 

2 2 + (y + _|_)2 


dr, + 2 (-4 f 

J z 2 + y 2 \ l 


dZ w (x ^ ) s eb (x) 

-r- d^ + 5 


For points on the wing surface, i.e. z = 0, s fc < y < s^: 


h d? 2 w (x,f J 


' (y + e, ) 2 (y + 

In 


V. K b 2 

^_)2( _ J2_) 

t l: 


+ /(-^ 


dR ]., i r £dz w (x >^ 


4 p i r 

dx y J 


. + V ^b 

J y ^ y 


"I d c 

V- 


. d 2 z w ( x,y ) r 

+ 2 d^ \ {s t 


+ (y - s t ) It, 


y s 


(Continued on next page) 



+ 


dZ (x,s ) ds 
W i 

dx dx 


gn |- 


(s„ - y) 2 (s # + y) 2 (y - -#7) 2 (y + ~s 


V 

S JJ 


V 2 


d Z,_,(x,s t ) ds t 


w 


dx dx 


in 


R 2 Rk 2 

(y- s )2(y+ s )2(y - -f-) 2 ( y + -f~) 2 

r _t S t I 

V 8 J 




+ 2 in (y 2 ) 


(50) 


v 


2 ,t 

U 7T 


i [ ft d V x >^> 

w lj 


dx l y + (- 


1 + — 1 + — -L 


y + 


V V 


d4 


y - 


/ dz w< x >£) dz v, (x >y) 


-/ 


s fl f — 2±— 
z \ dx 


dx 


Z - y 


az w u,y) - y 

at, 77 in — 


dx 


y - S . 


1 

+ 

y 




i dZ(x,£) S' (x) 

w ^ eb 

— d^ — + — r - 


: si ) 


w 2,t dz w (x,y) 

U dx 

00 



1 (52) 


We note that in this case no singularities occur on the body 
surface or in the gap between the wing and body. At the leading and 
trailing edges of the wing, however, the characteristic logarithmic 
singularity associated with two-dimensional incompressible flow at a 
sharp edge appears. 

As before, the thickness potential for the equivalent body cross 
section W 2 B is given by equation (17). "Because of the symmetry of 
these configurations, nonlifting flows will produce no lateral forces or 
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moments but only a longitudinal drag force. Since, as is the usual case 
of realistic wing— body combinations, the configurations considered herein 
have the body base located aft of the trailing edge of the wing tip 
chord (i.e. > X smQ ), the drag coefficient at zero lift is given by 
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eb 
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where Cp ^ -*- s t: ^ le ^ ra 9 coefficient of the equivalent body of revolution 

and is given by equation (22) . 

Elliptic bodies .- For the wing-body combinations with indented 
elliptic cross section and sweptback trailing edges, equation (36) for 
w 2 , t applies for x < X rte:L . For X rt6l < x < X snl2 , use of the 
Joukowski transformation, equation (30), provides, the result that 
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and again restrict 
equations (38) and 
components. For a 


attention to the first quadrant of the crossflow plane, 
(39) provide the following results for the velocity 
point at general location but not on the wing surface, 
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For points on the wing surface, i.e., z = 0, s fc < y < s ^ (or equivalently 

Zl = o, s lt < y x < s l£ ). 
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RESULTS AND DISCUSSION 


While experimental verification of the theory is considered essential, 
particularly at this stage of development, there is not available at the 
present time experimental data for transonic flows about wing-indented 
body combinations suitable for comparison with theory developed here. 
Although a parallel experimental program was originally considered for a 
typical member of the class of configurations described herein, that 
program has, unfortunately, been delayed. Cor sequently , experimental 
verification of the theory, particularly for pressure distribution compari- 
sons which are vital in assessing the validity of the assumptions of the 
theory within the various regions (body surface, wing surface, wing-body 
junction, wing leading and trailing edges, etc.) of the near flow field 
of these configurations will have to be deferred. 


In order to illustrate the general behavior of the theoretical 
results for transonic flows about the slender wing-body combinations 
considered here, the surface and flow field pressure distributions for 
several typical members of the classes of configurations described 
previously are given in figures 4 through 7. For example, in figure 4 
pressure distributions are presented for a finite thickness wing-indented 
circular body combination with a nonzero taper ratio and straight trailing 
edge in which the equivalent body is a parabolic-arc of thickness ratio 
D /l = 0.1, the wing is a truncated delta wing with aspect ratio AR = 1.7, 
and taper ratio TR =0.2 (so that 3^ = 58°, 3 = 0) , and has parabolic- 

arc profiles of thickness/chord ratio t/c = 0-04. The wing root chord 
is half the body length, i.e., C R // = 0.5, the root chord leading edge 


f = 0.86. The 


is located at = 0.25, and the body base is at 

longitudinal pressure distributions given in figure 4 are for the free- 
stream conditions M = 1, a = 0° and are presented at the two angular 
positions 0 = 0°, 90° in the crossflow plane and at locations on the 
body surface and also along lines parallel to the body axis but removed 
laterally from it by distances of 1, 2, and 4 times the maximum equivalent 
body diameter D. Thus, the pressure distributions given for O' = 0° 
and r/D =1, 2 cut across the wing surface, intersecting the leading edge 

at the axial positions x / P = 0.410 and 0.570, respectively. 
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The wing-body surface pressure distributions shown in the first plot 
of figure 4, when compared to the pressure distribution on the equivalent 
body alone, demonstrate the large effect that the wing has upon the body 
pressure distribution. Moreover, it clearly shows the rapid variation of 
the pressure distributions caused by the singularities at the points 
(x// = 0.320, 0.75) where the wing leading and trailing edges pierce the 
body surface and also at x/ / = 0.65 where the leading edge of the wing 
tip chord is. The discontinuities at the points where the leading and 
trailing edges intersect the body surface are related to the character- 
istic logarithmic singularity associated with the two-dimensional thick- 
ness problem (i.e. , <p 9 , ) of flow at a sharp edge. The discontinuity at 
the leading edge of the wing tip chord is due to the discontinuity in 
slope of the indented body at that point. This discontinuity occurs 

since in order that the equivalent body area distribution and its deriva- 
tives remain smooth, it is necessary for the indented body to have a 

slope discontinuity at x// = 0.65 to compensate for the one due to the 

wing. This discontinuity also occurs at the trailing edge of the wing tip 
chord and would be evident if the trailing edge were swept forward, so 
that the point (x = X ) would be separate and distinct from the axial 
location where the trailing edge pierces the body surface (x = X rtSi ) . 

For the case of a straight trailing edge as in figure 4, those points, of 
course, coincide so that only one singularity is evident. The flow field 
distributions shown at r/D = 1, 2, and 4 illustrate several interesting 
features. The most prominent is the propagation into the flow field of 
the singularities which occur in the surface pressure distribution at the 
three points discussed above. This is a direct consequence of using the 
transonic equivalence rule to provide flow-field information based upon 
knowledge of flow properties at the body surface. Also evident in the 
distributions along the lines r/D = 1, 2, 0 = 0° are the logarithmic 
singularities at x/.P = 0.410 and 0.570, respectively, as those lines 
cross the wing leading edge. . The longitudinal flow field pressure distri- 
butions provide insight into the rapidity with which the flow field 
becomes ax isymmetr ic and equal to that about the equivalent body. At the 
lateral distance r/D = 4, the pressure distributions at A = 0°, 90° 
at points ahead of the leading edge of the wing tip chord (X sm /£ = 0.65) 
are virtually indistinguishable from that about about the equivalent body - 
except for the exponentially small region of influence of the logarithmic 
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singularity propagating out from the body surface at the point where the 

leading edge pierces the body surface. However, within the axial region 

corresponding to the wing tip chord, 0.65 < x/ f < 0.75, the pressure 

distribution still shows some effect of the wing, although it is clearly 

diminishing. This is not surprising and could have been anticipated from 

the results from the delta wing with zero taper ratio given in figure 4 of 

reference 7 which also indicate that the effect of the wing on the flow 

field at lateral distances of several maximum body radii is negligible at 

all axial locations except those in the near vicinity of the wing maximum 

span. Knowledge of the region in which the flow about geometrically 

complex configurations of this type can be considered axisymmetric and 

equal to that about the equivalent body is quite significant and can 

provide, for example, useful information for a completely numerical 

finite difference solution in applying the far-field boundary condition. 

The drag coefficient for this configuration, which is provided by 

evaluating numerically the integral in equation (22) , is found to be 

C = 0.1044. 

Dt 

Analogous results are given in figure 5 for a lifting flow about 
this same configuration for the free-stream conditions M = 1 and 
a = 2°. We note again that the singularities discussed with regard to 
the nonlifting case also appear here. Moreover, due to the nature of 
lifting flows near a sharp edge, the logarithmic singularities associated 
with the thickness problem are further reinforced by the inverse square 
root behavior associated with the two-dimensional lifting problem 
(i.e., ) of flow around a sharp edge. The net result is the more 

rapid variation of pressure evident in those regions. Nevertheless, the 
flow field distributions again display the strong tendency to return to 
those generated by the equivalent body alone, as is most apparent in the 
flow field distributions at r/D = 4. At this angle of attack, equations 
(27) , (28) , and (29) provide the following results for the aerodynamic 

coefficients: 

C = 1.7070, C = 0.1342, C = -0.8906 
1j u ’ m 

In order to demonstrate the pressure distribution behavior typical 
of the wing-body combinations considered here having swept-back trailing 
edges and non-zero taper ratios, results are given in figure 6 for a finite 
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thickness wing-indented circular body combination in which the equivalent 
body is a parabolic-arc of thickness ratio D // = 0.10, the wings have 
parabolic-arc profiles of thickness ratio t/c = 0.04, planform aspect 
ratio AR = 2.8, taper ratio TR = 0.4, root chord C R ^/ i = 0.3, with the 
root chord leading edge at = 0.25 (so that d^ e - 45°, 6 fce = 23.75°). 

Analogous results are presented in figure 7 for a finite thickness wing- 
indented elliptic body combination composed of a parabolic-arc body of 
semimajor to semiminor axes A = 3 and a wing essentially identical to 
the one described above for the circular body except that the trailing 
edge is swept at the angle 6^ = 25.05°. The trailing edge sweep angles 

of these configurations are such that the axial locations of leading edge 
of the wing tip chord and the point where the trailing edge pierces the 
body surface coincide, i.e., X smi = X rt6i . 


In figure 6, we note that the general variation of both surface and 

flow field pressure distributions are essentially the same as those of the 

straight trailing edge configuration shown in figure 4 for points ahead 

of the leading edge of the wing tip chord (X sm ^/I = 0.571). However, 

within the axial region containing the wing tip chord < x < X am _) 

o rn b m 2 

and coincidentally, the wing trailing edge, the pressure distributions 

now indicate a much more rapid variation, while still exhibiting the same 

trend as that shown in figure 4. Apparently the gap between wing and 

body in the crossflow plane within this region influences the behavior of 

the surface pressures to a greater degree than in the case when there is 

no gap, and consequently an unbroken lateral distribution of sources 

along the wing from body to wing tip. The flow field distributions within 

the near flow field of this region maintain this rapid variation, with the 

distributions of C at 9 = 0°, r/D = 1, 2 also exhibiting the character- 

P 

istic influence of the logarithmic singularities at the points where 
these lines cross the trailing edge. Nevertheless, beyond the wing tip 
the flow still displays the characteristic tendency to return to that of 
the axisymmetric flow about the equivalent body, as is evident in the 
distribution at r/D = 4, a distance which is only slightly beyond the 
point of maximum span, r/D = 3.2. 


The surface and flow field pressure distributions shown in figure 7 
for the wing-indented elliptic body combination described above are 
essentially similar in behavior to those in figure 6 for the corresponding 


36 



circular body and do not exhibit any new characteristic features. We note 
that the asymmetry introduced by ellipticity of the cross section alone 
(excluding the influence of the wing) , while being evident at the body 
surface, rapidly dies out. In reference 6, it was shown that for a smooth 
elliptic body alone having a semimajor to semiminor axis ratio A = 3 the 
flow field becomes essentially axisymmetric at r/D = 1. 

Perhaps the most notable feature of the theoretical results presented 
here (and in ref. 7) for transonic flows about the classes of wing-indented 
body combinations being considered is the behavior of the pressure distri- 
butions caused by the singularities which occur at the following axial 
locations: 

• X . — leading edge pierces body surface 

r '- e x 

• — trailing edge pierces body surface 

• X sm — leading edge of wing tip chord 

• X gm — trailing edge of wing tip chord 

Although the singularities at general locations along the leading and 
trailing edges, for either nonlifting or lifting situations, could be 
included here, they are not, both since their character and origin are 
well known and also because they are local phenomena and, consequently, 
of restricted influence unlike the singularities delineated above. 

It is important to realize that the basis of these singularities is 
essentially geometric in character, with the difficulty arising from 
either a discontinuity in first (at x = X sm ^ , X gm ) or second 
(at x = X r ^ e , X r £ e ) derivative of the indented body area distribution, 
which causes, in turn, discontinuities in the surface velocity components. 
Then, because the transonic equivalence rule is used to provide flow field 
information based upon knowledge of flow properties on the body surface, 
these discontinuities are propagated laterally into the flow field making 
their presence even more evident. A direct method of alleviating this 
problem, while at the same time providing both a more general and 
realistic approximation would be to smooth these junction points with 
monotonically varying fillets. It appears, however, that a simple func- 
tional representation of fairing curves of this nature is not possible. 
Analytic (i.e., cubic), trigonometric, or exponential curves, while 
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satisfying the end conditions of matching slope and ordinate at two points, 
fail to be continuously monotonic under boundary conditions typical of 
the configurations considered here. Thus, a wiggle would result in the 
faired curve and this is unacceptable. A means of eliminating this 
problem would be with piecewise continuous spline-fit functions. In any 
case, by whatever means the smoothing is accompl ished , the result would 
be a more accurate representation of the actual solution in the vicinity 
of these points. 


CONCLUDING REMARKS 

Theoretical analysis and development of associated computer programs 
have been conducted in order to develop calculative techniques for pre- 
dicting properties of transonic flows about certain classes of slender 
wing-body combinations. The theoretical analysis is based upon a 
combination of the transonic equivalence rule and uses either an arbi- 
trarily specified solution or the local linearization method for deter- 
mining the nonlifting transonic flow about the equivalent body. 

Computational programs, which are documented in a general user’s 
manual and included as part of this report, have been developed for 
finite thickness wing-body combinations in which the bodies are area- 
rule indented in such a manner that the resultant equivalent bodies 
remain smooth. The equivalent body profiles are either user-supplied 
subject to certain continuity and closure restrictions or program-supplied 
in which case the radius is of the general class R ~ x / l - (x/f) n or 
1 - x/ 1 - (1 - x//) n . In addition, the body cross sectional shapes are 
either (1) circular or (2) elliptic and such that a constant ratio A of 
semimajor to semiminor axes is maintained along the entire body length. 

A general class of wings is considered which are symmetric in plan- 

form about the azimuthal body meridian (x-z plane) and consist of 

straight leading and trailing edges swept at arbitrary angles. The 

positions of the leading and trailing edges of the root chord are located 

at arbitrary locations on the body axis, and the profiles are described 

by ~ x/c w - (x/c^)™ or 1 - x/c^ - ( 1 - x/c^) m where x is the 

axial distance from the leading edge and c is the local chord. 

w 
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These programs provide longitudinal pressure distributions for both 
nonlifting and lifting situations, at arbitrary angular positions in the 
crossflow plane at points along the body and wing surface and also along 
lines parallel to the body axis but removed at arbitrarily specified 
lateral distances from it. In addition to the pressure distributions, the 
aerodynamic characteristics of lift, drag, and pitching moment are also 
provided. 

The theoretical pressure distributions predicted by these programs 
for certain members of the class of configurations described above indi- 
cate quantitatively the relatively large effects of wing thickness and 
lift on both the body and flow field pressures, and also serve to point 
out the singularities inherent to the theory as it is presently constituted. 
In addition, they demonstrate the large influence that sweeping the trailing 
edge and introducing a finite tip chord has upon the pressure distributions. 

In conclusion, we emphasize that the techniques employed here are 
quite fundamental and possess great generality so as to allow extension 
to even more complex configurations. Moreover, since the solutions to the 
various two-dimensional crossflow problems are independent of Mach number, 
they can be calculated once and for all once the geometry of the configura- 
tion is fixed and then combined with any one of a possible variety of 
solutions (experimental, numerical, etc.) for the transonic flow about the 
nonlifting equivalent body. We suggest, furthermore, that experimental 
work be conducted to determine surface and flow-field pressure distribu- 
tions on selected wing-body combinations in order to define more clearly 
the extent to which the theory applies to configurations of this nature 
and, also, that consideration be given'to developing methods to smooth 
the solutions in the vicinity of the various discontinuities in the area 
derivatives of the indented body shapes. 


Nielsen Engineering 
Mountain View, 
July 13, 


& Research, 

California 

1972 


Inc. 
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APPENDIX A 


COMPUTER PROGRAM USER’S MANUAL 
SUMMARY 


An operating manual is given for the computer program developed in 
conjunction with the theoretical work presented in this report. The 
program computes the transonic surface and flow field pressure distri- 
butions and aerodynamic characteristics for various classes of wing-body 
combinations considered herein. Use is made of the transonic equivalence 
rule and either the local linearization method or a user-supplied 
solution for flow about the nonlifting equivalent body. 

A description of the general operating procedure of the program is 
given, together with instructions for the preparation of input data, 
sample output of test cases, and a listing. The program is written in 
FORTRAM IV programming language and prepared specifically for use on an 
IBM 360/67 series computer. Typical running times are approximately 
30 to 45 seconds for the equivalent body calculations using the local 
linearization method and about 2 minutes for the crossflow solution 
calculations involving approximately 1000 points located typically 
along the wing, wing-body junction, and flow field. 

DESCRIPTION OF PROGRAM 

The computer program presented here is applicable to several 
classes of finite thickness wing-body combinations discussed in the 
preceding report in which the bodies are area-rule indented along the 
wing-body junction in such a manner that the total cross-sectional area 
distribution is identical to that of a smooth body having a specified 
profile. The programs compute the surface and flow field pressure 
distributions, for both nonlifting and lifting situations for straight or 
swept forward trailing edge planforms and for nonlifting situations for 
swept back trailing edge planforms, at arbitrarily specified angular 
positions in the crossflow plane, at points along the body and wing sur- 
face, and also along lines parallel to the body axis but removed from it 
at specified lateral distances. In addition, the aerodynamic character- 
istics of lift, drag, and pitching moment are computed. 
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V. 

APPENDIX A 


Wing and Body Geometry 

The wing and body geometries of the configurations programmed are 
shown schematically in figures 2 and 3. Figure 2 illustrates two 
members of the class of wing-body combinations which have indented 
bodies that are circular in cross section, while figure 3 shows the 
corresponding members of the class having indented bodies with elliptic 
cross section. 


Program-supplied equivalent body profiles . - Unless the user specifics 
to the contrary, the class of equivalent bodies of revolution of both 
types of the above configurations consist of profiles described by the 
equations 
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with n = constant 2 . In reference 7, the profiles of the elliptic 
bodies were restricted to parabolic arcs, i.e. equations (64) or (65) 
with n = 2. Thus, this work extends the elliptic body category to 
include the entire class of equivalent body profiles used for the case 
of the circular bodies. 

User-supplied equivalent body profiles .- At the user's option, an 
arbitrarily specified equivalent body profile may be substituted in 
lieu of the above class of profiles. The modifications necessary to 
the program are deta’iled in the PROGRAM INPUT section. The restrictions 
on these profiles depend, in part, on the method used to calculate the 
solution for the nonlifting flow about the equivalent body, i.e. u . 

If the local linearization method, as presently constituted, is used 
to calculate u 0 (see eqs. (5), (6), and (7)), then it is necessary 

that the profiles be closed, have sharp tips, and have continuous 
derivatives through the fourth. On the other hand, if the solution for 
.Uj. is user-supplied, then the requirement from the other portions of 
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the solution, i.e. 02, a > 02, t> anc3 02, B> -*- s the equivalent body 

profiles have continuous derivatives through the second. 

Indented-body profiles .- The ordinates of the indented body profiles 
are fixed once the equivalent body profile and wing profile are specified. 
For circular bodies with straight/swept forward trailing edge planforms, 
the indented body radius R^ is found through a Newton-Raphson iteration 
procedure on the expression 


* R l h = * R 2 + 4 J Z w (x,fJ d£ 

«b 


( 66 ) 


while the derivatives dR^/dx and d 2 K^/dx 2 are calculated by using an 
appropriate five-point difference formula. For sweptback trailing edge 
planforms, the above method applies up to x = X r te L • For x rtei < x < x sm 2 » 
R^ is found without iteration from the expression 


IT 



ir R 2 + 4 



V x »*> 


di* 


(67) 


Analogously, for elliptic bodies with straight/sweptforward trailing 
edges the semimajor axis a of the indented elliptic cross section is 
found by iteration on 
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R 2 

eb 
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s 

J z w (x,4)'d$ 
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( 68 ) 


with the derivatives da/dx and d 2 a/dx being evaluated numerically by 
using the appropriate five-point difference formula. Fow sweptback 
trailing edge planforms with X rte < x < X sm2 , a is found directly 
from the expression 


Tra 
eb X 


tt R 2 = 




Z w (x,^) d| 


(69) 
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The general class of wings considered for both types of body 
shapes described above have wing planforms that consist of symmetric 
straight leading and trailing edges, swept at arbitrary angles i3 p and 
(3 te , respectively, to the y axis. Both p3^ g and 3 te > are measured 
positive clockwise; thus, for (3 less than, equal to, or greater than 
zero, the trailing edge is correspondingly sweptforward , straight, or 
sweptback. The position of the leading edge of the wing root chord 
X rfe anc ^ -*- ts •'- en 9 ttl c Rt are ar bitrary. The wing profiles are 
represented by expressions of the form 


Z 

w 

c 

w 


(m/m-i) ✓ 

w ( x 

2( m - !) [c^ 



(70) 


or 

(■ - ©") 

where is the local chord, x the distance from the leading edge, 

m is a constant 2 . 2 , and t w is the wing thickness-to-chord ratio. 

The wings are assumed to maintain a constant thickness-to-chord ratio 
across the span, with the consequence that the wing profiles at all 
spanwise locations are geometrically similar so that 


(m/m-i ] 


Z t m 
w _w 

c ~~ 2 (m - 1 ) 

w ' 



\l = (2 w (x >y )} 

2 ' c w<V> 


(Z (x , o) ) 
w * max 


'R. 


(72) 


General 

The coordinate system used in the program is a body-fixed Cartesian 
system centered at the body nose with the x axis directed rearward and 
aligned with the longitudinal axis of the body, the y axis directed 
to the right facing forward, and the z axis directed vertically upward, 
as shown in figure 1. Because the transonic equivalence rule allows the 
perturbation potential 0 to be expressed in the form 
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0 ~ 02 , ct + 02 , t: ~ 02 , B + 0p (7 3) 

where each of the components has the meaning indicated in figure 1 , 
and since 0 2 q, 0 2 and 0 2 3 satisfy the two-dimensional Laplace 
equation 

(02, i) + (0 2 , i) = 0 . (74) 

yy zz 

they are independent of Mach number. Consequently, once the geometry of 
the configuration is fixed they can be calculated once and for all, 
stored, and used, for example, in a comparative study of a certain wing- 
body combination as the Mach number is varied systematically throughout 
the transonic range. An option for doing this is available and is 
discussed in the PROGRAM INPUT section. The only portion of the solution 
dependent upon M^ is 0 ^ and this term represents the solution to the 
full transonic equation ( 1 ) for flow about the nonlifting equivalent 
body. 

Local linearization solution of u .- If the local linearization 
method is used to determine the solution for 0 , or more conveniently, 

U B = (0q^ x ’ t ^' en accor ding to whether M^ is near one, below the lower 
critical, or above the upper critical, equation (5) , ( 6 ) , or (7) must 

be integrated. Since these are all first order ordinary nonlinear 
differential equations, appropriate initial conditions are required. 

These are given at the point x , which is the positive root of the 
equation 

S^ b (x) = 0 (75) 

that is closest to the origin. The values of Ug/U^ at this point 
are, for accelerating transonic flows with M^ ~ 1 (eq. (5)) 


M 2 (y + 1) 

OQ I 



X 



S" (x) 
oh 

' X 





for purely subsonic flow (eq. ( 6 )) 


( 76 ) 
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U B _ I f 

U 4t t I 

00 -i 


s ; b (x) 


s ; b ( ^ 


x - ^ 


d£ 


and for purely subsonic flow (eq. (7)) 



o 


(77) 


(78) 


The integrations start at x , proceed to a specified point near 
the nose, and upon reaching that point, return to x^, restart the inte- 
gration procedure, and then continue toward the tail. In each of these 
programs, the differential equations are integrated b y using Hamming's 
modified predictor-corrector method described in the Scientific 
Subroutine Package (SSP) available from the IBM Corporation. The 
integrals involved in those differential equations are evaluated by 
using Simpson's rule. 

User-supplied solution for u n .- At the user's option, an arbitrarily 

B 

specified solution for u n can be used in lieu of the local linearization 
solution. This solution can involve mixed transonic flows with imbedded 
shocks and can be determined in any of a variety of ways (numerical, 
experimental, etc.). Details regarding the manner of inputting this 
information to the program are discussed in the PROGRAM INPUT section. 

Crossflow Potentials and Aerodynamic Characteristics 

This section assembles for user convenience, the crossflow potentials 
and aerodynamic characteristics of all of the configurations considered 
in this report. 

Straight/Sweptforward trail ing-edge planforms .- For the classes 
of finite thickness wing-circular body combinations considered herein 
which have straight/sweptforward trailing edge planforms, (see fig. 2(a)), 
the following results are provided for W 2 a , W 2 ,t> an< 3 W 2 ,B at the 
indicated axial locations: 
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U 2ir 

OO 


In a (0 < x < £.) (83) 


where s in equation (82) denotes either the leading (s = s^) or 
trailing (s = s^) edge, depending upon the axial location. We note that 
for all of the configurations considered in this report, the solution 
for W 2 g is given by equation (83). 

The aerodynamic characteristics of this class of configurations are 
given as follows. The drag coefficient at zero lift is found through 
a numerical integration of the expression 
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dS eb (x) 
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( 84) 


where C D ^ as the dra< 3 coefficient of the equivalent body alone, S m 
is the maximum area of the equivalent body, and Cp ^ is the pressure 
coefficient on the surface of the nonlifting equivalent body and is 
equal to 


__B / /dR eb (x) 

U oo \ dX 


(85) 


Because of the symmetry of these configurations no lateral forces or 
moments exist at a = 0. For the lifting situation, the coefficients 
of lift, drag, and pitching moment are given by 



( 86 ) 



a=o 
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where X srn2 and arG the axial locations, respectively, of the 

trailing edge of the wing tip-chord and the point where the wing leading 
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edge pierces the body surface. The integrals involved in evaluating 
the pitching moment are calculated by using Simpson's rule. 

The corresponding results for the wing-elliptic body combinations 
(see fig. 3(a)) are for W 2ja , W 2>t : 
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2 ,a 
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X rie 
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where s x in equation (92) denotes either the leading (s-l = s i^) or 
trailing (Sj^ = s lt ) edge in the transformed o x plane. 

The aerodynamic characteristics of these elliptic body configurations 
are given, for nonlifting flows, by 
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- ir In R 



( 93 ) 
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where c D e }-> is given by equation (84) and K(. : ,) is the complete 
elliptic integral of the first kind; and for lifting flows by 



(96) 


where the integrals involved in evaluating the pitching moment are 
calculated by using Simpson's rule. 

Sweptback trailing edge planforms .- For the classes of finite 
thickness wing-circular body combinations having sweptback trailing 
edge planforms considered here (see fig. 2(b)), the following results 
for w 2 , t are provided: 
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and the drag coefficient at zero lift is given by 
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a=o eb 


( 100 ) 


where C D n is given by equation (84) 
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The corresponding results for the wing-elliptic body combinations are 
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with the drag coefficient at zero lift c D a=0 given by equation (93) . 


Operating Procedure 

The basic characteristics and general operating procedure of the 
computer program developed herein are straightforward and can be out- 
lined as follows. After reading the input data (which is detailed in 
a subsequent section) and checking it for obvious errors, the program 
proceeds to calculate certain required geometrical and flow-field 
constants. Then, if the user selects the equivalent body profile to 
be of| the class described by equations (64) or (65), the program proceeds 
to calculate the exponent n from information regarding the point of 
maximum thickness (see equations (12), (14)). The point x g is next 

found by solving equation (75). If however, the equivalent body profile 
is user supplied, the calculation of n and x^ are omitted. Next, 
the axial locations X r |j e ^ and X r ^- ei , which represent, respectively, 
the points where the wing leading and trailing edges pierce the body 
surface, are calculated. The exponent m describing the wing ordinates 
(see eqs. (70), (71)) is then calculated in a manner similar to that 
used to determine n. The calculation of n, x g , X r ^ Si , X rtei , and m 
are all performed in an iterative fashion by using the standard Newton- 
Raphson iteration scheme. 

With the calculation of the above parameters complete, the program 
prints a number of geometrical and flow-field characteristics for the 
case at hand. If the solution of u Q is user-supplied, the program 
begins at a point close to the nose (x/f = 0.005) and proceeds toward 
the tail. If the local linearization method is used to determine u n , 

O 

then the appropriate initial value (eqs. (76), (77), or (78)) for the 

local linearization equation at hand (eqs. (5), (6), or (7)) is 
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calculated at x = x g and the numerical integration begun. In the case of 

purely sub- or supersonic (eqs. (6), (7)) flow, it is convenient to redefine 

the dependent variables (see eqs. (80) through (88), ref. 6) and integrate 

a simplified differential equation. For the M w ~ 1 case, it is more 

advantageous to integrate equation (5) as it stands. Because of the 

special character of that equation at x = x , it is necessary to use 

a Taylor series for u^/U^ in the neighborhood of that point in order 

to avoid a singularity in the numerical integration. Consequently, for 

that case in addition to u„/U several derivatives are also required 

B 00 

and are calculated by the program. Details are given in reference 6. 

The numerical integrations then continue toward the nose and stop at 
a point (x/f = 0.005) close to it. The integrations are not carried 
directly to the nose because, although this is possible for the purely 
supersonic case, the local linearization method predicts a logarithmic 
singularity at x = 0 for a sharp-tipped body, much like that indicated 
by linearized theory. With the integration to the nose complete, the 
program returns to x , restarts the numerical integration, and continues 
toward the tail. As these calculations progress (using either a user- 
supplied or local linearization solution for u ) , the surface and flow- 

■D 

field pressure distributions are calculated from equations (2) and (3) 
and the output printed at specified axial locations. Until the point 
x = X r ^ e ^ is reached, the appropriate crossflow solutions for determining 
these pressure distributions are those for the smooth body alone (see 
eqs. (79), (81), (83), (89), (91), and (101)). Beyond X r ^ Sl , for the 

case of nonlifting flows, the crossflow solutions are calculated from 
equation (82), (83), or (92) and (83) for x r ^ 6l < x < X rtei . Beyond 

X r te x > f° r planforms with straight/ sweptforward trailing edges, the 
crossflow solutions revert to those for the smooth body alone; while 
for planforms with sweptback trailing edges, the appropriate crossflow 
solutions for X r t 6l < x < X sm2 are given by equations (99) , (83) , 
or (103) , (83) and beyond X sms the solutions revert to those for the 
smooth body alone. The calculations continue until the body base is 
reached, i.e. x = X^; the calculation then returns to the mainline pro- 
gram, prints the value of the drag coefficient, and reads the input 
data for the next case. For lifting flows about planforms with 
straight/ sweptforward trailing edges, the calculations proceed in a 
similar fashion to that of the nonlifting case (with the additional 
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output of surface and flow-field pressure distributions at the angular 
locations i 9 rather than just + 9 ) until the axial location of the 
trailing edge of the wing tip chord is reached, i.e. x = X sm ^ . Beyond 
that point, for reasons given in reference 7, no further pressure distri- 
butions are given. However, the calculation of flow about the equivalent 
body, i.e. u Q , continues to the body base in order that the drag 
coefficient at zero lift c D a=0 can determined. When x = X^, the 
calculation returns to the mainline program, determines the coefficients 
of lift, drag, and pitching moment from equation (86), (87), and (88), 
or (94), (95), and (96), prints these values, and then proceeds to 
read the input data for the next case. 

PROGRAM INPUT 

The variables that are input to the program are described in the 
following list; 


AL 

ALPHA 

AMACH 

ANGLE 

CRT 

ICOPY 


MAREA 


MOPT 


Dictionary of Input Variables 

ratio of semimajor to semiminor axis (a/b) of elliptic 
cross section: program default value AL = 1 

angle of attack, in degrees; program default value, 
ALPHA =0 

free-stream Mach number 

sweep angle, in degrees, of wing leading edge (measured 
positive clockwise from y axis and restricted to 
values 0 < ANGLE < 90, see figures 2, 3) 

wing root chord normalized by complete body length, 


integer index for program option for using previously- 
stored values of crossflow solutions; equal to 0 or 1; 
program default value, ICOPY = 0 

integer index for program option for using user-supplied 
or program-supplied subroutines for equivalent body area 
and derivatives, equal to 0 or 1 ; program default value, 

MAREA = 0 

integer index for program option for using user-supplied 
distribution for u g or program-supplied local linearization 

solution; equal to 1, 2, 3, or 4. 
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NTH ETA 
NXEB 

RF ( I ) 

SSMAX 

TAUB 
TAUW 
THETA (I) 

TR 

XEB (I) 

XLBASE 

XLOUTP 

XRLE 

XMTB 

XMTW 

XS2EB 


integer indicating number of angular positions in crossflow 
plane at which output is desired; 1 ^ NTHETA 5 


integer indicating number of table entries of the user- 
supplied distribution of ; NXEB 201 


six-dimensional vector representing values of r/D (the 
radial distance in the crossflow plane normalized by the 
maximum equivalent body diameter D) at which flow-field 
pressure distributions are to be calculated 

maximum wing semispan normalized by complete body length, 

s /[ 
max 


thickness ratio of equivalent body, D / l 


thickness-to-chord ratio of wing (see eq. (10) ) 

NTHETA-dimensi-onal vector representing values of the angle 
0 (in degrees) in the crossflow plane; 0 <_ THETA (I) 90 


wing planform taper ratio, 


■WV 


0 < TR <r 1 


NXEB-dimens ional vector representing values of the axial 
locations (normalized by the complete body length) where 
values for the user-supplied distribution u fi are given 

axial location of body base normalized by the complete body 
length, X^/f 


interval size, as fraction of complete body length, 
between output stations for pressure distribution print- 
out 

axial location of leading edge of wing root chord normalized 
by complete body length, X^^/f. 

axial location of position of maximum thickness of equivalent 
body of revolution normalized by complete body length 
(see eqs . ( 12) , ( 14) ) 

location, as fraction of distance from wing leading edge 
to local chord length (x/c^) , of position of wing maximum 

thickness (see eqs. (8), (9)) 

user-supplied value of the axial location, normalized by 
the complete body length, where the user-supplied equivalent 
body profile satisfies S” b (x) = 0; only necessary as 

input when user supplies equivalent body profile and also 

uses local linearization method to calculate u D 

B 
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UEB(I) NXEB-dimensional vector representing values of Ug/u^ for 

the user-supplied distribution of u D 

.D • 


Input Format and Options 


All of the input variables are entered into the program under a 

NAMELIST format (the one exception being the vectors UEB(I), XEB(I) which 

represent the ordinates and abscissas, respectively, of the user-supplied 

velocity distribution u_ and the format of these quantities is dis- 

B 

cussed below) . The name of the NAMELIST data block is TRANIN. 

Default values .- The following input variables have default values 
that are indicated below and unless the user wishes to change them, it 
is not necessary to enter them in the input data block 


Variable Name 


AL 

ALPHA 

RF(I) , 
1 = 1 , 2 , 3 , 
4 , 5,6 

NTHETA 

THETA ( I ) , 

1 = 1,2 


Default Value 

1 . 

0 . 

1 ., 2 . , 3 . , 4 . , 5 . , 6 ., 


2 


0 . , 90 . 


MAREA 

ICOPY 

XLOUTP 


0 

0 

01 


It is important to realize that the above variables assume their 
default values each time the program is run. If the user wishes any 
of the above variables to be different from its default value, this 
must be specified in the data statement for each run. All other input 
variables, once specified, remain unchanged by the program; thus, it is 
unnecessary to respecify them in subsequent runs if their values are to 
remain constant. 

Local lineariztion option .- To use the local linearization method 
to determine u Q , it is necessary to specify in the input data the 
appropriate value for the integer index MOPT. Depending on the free 
stream Mach number, the proper value of MOPT to use the local lineari- 
zation method is 
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Free Stream Mach No. MOPT 

~ 1 (near sonic flow) 1 

M < M „ (below lower critical) 2 

00 — cr , £ 

M < M (above upper critical) 3 

cr.u — 00 


User-supplied u_ option.- If the user wishes to supply the 
solution for Ug, then the program will bypass the local linearization 
calculations by specifying MOPT = 4. The solution for u B is read 
into the program in the form of a tabular input of values of 
u g/ U oo “ vs - x/i immediately after the NAMELIST input block. Provision 
has been made for inputting ordinate and abscissa values up to a 
maximum number of 201 each (UEB(201), XEB(201)), i.e. values at each 
half percent of the body length if equally spaced. It is assumed that 
a sufficient number entries are made that linear interpolation in the 
table is appropriate. 

User-supplied equivalent body profile .- If a class of equivalent 
body profiles not included in equation (64) or (65) is desired, then 
the user must set the integer index MAREA = 1, remove the following 
function subroutines from the program, 

•FUNCTION SEBPI(DZ) 

•FUNCTION SlEBPI(DZ) 

•FUNCTION S2EBPI (DZ) 

•FUNCTION S3EBPI (DZ) 

• FUNCTION S4EBPI ( DZ ) 

and replace them with his own. The above subroutines which are non- 
dimensionalized by normalizing them with respect to the body length, 
are defined in the following fashion. 

S eb (x ^> 

— 4 ~[2 = SEBPI (x/ i ) 


4tTP 


1 dS (x/i) 

¥F— /TP = slEBPI(x/i) 
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s ; b (x/i) 

4ttP 


d2S eb (x/i) 

4~| 7 d ('x7'i') g " S2EBPI (x/£) 


S^(x/i) 

4irf 2 


1 d3s e -b (x/i) 
47Ti 2 d (x/f ) 3 


= S3EBPI (x/ i) 


S 'V (x/l) 
4w£ 3 


1 d4s eb {x/i) 

47F ~ (x/|,r = S4EBPKX/I). 


We note again the requirements that if the local linearization method 
is to be used with these subroutines, then the functions must be such 
that the complete profiles are closed, sharp-tipped, and have continuous 
area derivatives through the fourth. In addition, the user must 
supply the point x/f = XS2EB, i.e. the point closest to the origin 
where 

S eb ( x / ^ ) * 0 

If, however, the user supplies both the equivalent body profiles and 
the distribution of Ug/U^ - vs - x/f, then it is only necessary that 
the equivalent body area derivatives be continuous through the second. 

Also, for this case, it is unnecessary to specify XS2EB. These require- 
ments are summarized below when user-specifying the equivalent body profile: 


Value of MOPT 


Requirements 

MOPT = 1,2,3 

© 

S S * -i , S". , S"; , S ’Y continuous 

eb * eb ’ eb * eb * eb 

for 0 < x < £ 


© 

User must input value of x/f =XS2EB 
where (XS2EB) = 0 


© 

Set index MAREA =1 in input 

MOPT- = 4 

© 

S , , S' , S" continuous for 

eb ’ eb ’ eb 

0 < x < £ 


© 

Set index MAREA =1 in input 
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An illustrative example of a user-supplied equivalent body 
profile and the corresponding area and derivative subroutines to be 
used with, say, the local linearization method (MOPT = 1, 2, or 3) can 
be given as follows. Consider an equivalent body profile formed by the 
top half of the sinusoidal curve given by 


r/J 


x / 1 



so that 


R eb (x/i) T eb . x. . 

J = -j- sin (v (j) ) 


with say t = .1 


S eb (x/?) T lb 




16 


(1 - cos (2 u (x /£) ) ) 


S eb<*/» 
4m f 2 


UT eb 


sin ( 2 tt ( x/ £ ) ) 


47 T (■ 2 


ir 2 t 2 , 
eb 


cos ( 2 v (x /£) ) 


S eb (x "> 
4ir £ 2 




sin ( 2ir (x/f ) ) 


S 'V (x/£) 


~ 7r ' t Tg- b COS ( 277 (x/f) ) 


The function subroutines for SKBPI(x/f) and SlEBPI (x/ £) with 
T = .1 are given by 

FUNCTION SEBPI(DZ) 

TAU = . 1 

PI = 3.1415027 

SEBPI = TAU*TAU* ( 1 . - COS (2 . *PI*DZ) )/16. 
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RETURN 

END 

FUNCTION SIEBPI(DZ) 

TAU = .1 

PI = 3.1415927 

SlEBPI = PI *TAU*TAU*SIN ( 2 . *PI *DZ ) /8 . 

RETURN 

END 

The FUNCTION subroutines S2EBPI, S3EBPI, and S4EBPI are given 
in analogous fashion. The point where 


S; b (x/J) = 0 


is given by 


XS2EB = 0.25 

and must be included in the NAMELIST input data block. 


Repetitive calculation storage option. - If the user wishes to 
undertake a systematic study of a wing-body configuration of the classes 
considered herein in which the geometry of the configuration is frozen 
and the Mach number and/or angle of attack are varied, a provision is 
included in the program whereby the velocity components associated with 
the crossflow solution for w 2 t - which is independent of M^ and 
a, and is by far the most time consuming crossflow potential to calculate 
is stored at the user-specified output locations for later use. These 
velocity components are then provided for the remainder of the cases to 
be run rather than recalculated unnecessarily. 

In order to activate the repetitive calculation option and make 
use of previously-stored results from a base run, it is necessary to 
set the integer index ICOPY = 1. The default value for ICOPY is 
ICOPY = 0 and this default value instructs the program to perform the 
crossflow calculations at the user-indicated axial locations and then 
automatically store these results in anticipation of use with the next 
case. If it is not desired to use those stored results for the next 
run, (i.e. ICOPY = 0 for the next case) the program simply replaces 
the previously-stored results with the ones being currently calculated. 
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Because for lifting situations, the crossflow calculations only 
proceed to x/i = X gm2 (as opposed to the a = 0 case, which carries 
the calculation to x / £ = X^) unless all of the cases in the study are 
for a / 0, the initial or base run which stores the crossflow results 
should be made for a = 0. 

Finally, when using the storage option, it is necessary, because of 
the different starting conditions involved, to use the same method for 
calculating u Q , i.e. (1) local linearization or (2) user-supplied 

distribution of u„ . 

B 


Data Format 

The data format is most easily demonstrated by an example. 

Consider the case of a wing-elliptic body combination composed of a 
parabolic-arc equivalent body of thickness ratio 1/10, the ratio of 
semimajor to semiminor axis of the elliptic cross section is 3, the 
body base is at 85 percent of the complete body length, the wing 
profiles are parabolic _ arcs having a thickness/chord ratio of 0.04, 
the wing root chord is 40 percent of the complete body length with the 
leading edge of the wing root chord located at x/£ = 0.3, the leading 
edge swept at 45 degrees, a taper ratio of 0.3, and a wing semispan being 
28 percent of the complete body length (this implies the trailing edge 
is straight) . The pressure distributions are required to be output at 
every 2 percent of the complete body length at angular locations 
9 = 0°, 45°, 90°, at radial distances r/D = 1, 1.5, 2, 2.5, 3, and 3.5 
in the crossflow plane, at sonic free-stream conditions and 2 degrees 
angle of attack, by using the local linearization method to 
determine the flow about the equivalent body. 

Thus, the input data cards would read (note that with a NAMELIST 
format, input variable sequencing is arbitrary) : 


CARD NO. 1 


COLUMN NO. 

2 9 

80 


& TRANIN 

AL= 3 . . AMACH= 1 . ,MOPT= 1 , TAUB= . 1 . TAUW= . 04 . XMTB= . 5 , 

CARD NO. 2 



COLUMN NO. 

2 

80 


XMTW= . 5 ,XRLE= . 3 ,CRT= . 4, ANGLE=45 . ,TR= . 3 , SSMAX= . 28 , NTHETA= 3 , 
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CARD NO. 3 


COLUMN NO. 

2 80 


THETA (1)=0. .THETA (2) =40. , THETA ( 3 ) =90 . , RF ( 1 ) = 1 . , RF ( 2 ) = 1 . 5 , 

CARD NO. 4 


COLUMN NO. 

2 80 


RF ( 3 ) = 2 . , RF ( 4 ) = 2 . 5 , RF ( 6 ) = 3 . , RF ( 6) =3 . 5 ,ALPHA=2 . .XLOUTP=.02. 

CARD NO. 5 


COLUMN NO. 

o 

00 

CM 


&END 


If for this case the user wished to supply his own distribution of 
Ug/Uoo - vs. - x/i, this would have been done by specifying MOPT = 4 
and also the integer NXEB representing the number (say, for example, 

101) of values of tu/U^ being entered (i.e. NXEB = 101) in the NAMELIST 
input data block. Next, all of the NXEB (101, in this example) values of 
u /U = UEB(I) would be read in under the card format 8F10.0, with each 
successive value of u_,/U occupying a space of 10 columns (including 
decimal point) with 8 values per card. Finally, the NXEB values (101, in 
this example) of the corresponding axial locations x / £ = XEB(I) of the 
above values of u^/U would be read in under the same format. Thus, 

Card format for UEB(I): Format (8F10.0), decimal point required 



Card format for XEB(I): Format (8F-10.0), decimal point required 


COLUMN NO. 
Data for 


10 

20 

70 

80 

XEB (1) 

XEB (2) 


XEB (8) 


I etc . 
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COLUMN NO. 
Data for 


10 

20 

30 

40 

50 ( 

XEB {97) 

XEB (98) 

XEB (99) 

XEB (100) 

XEB (101) J 


MESSAGES PRINTED BY THE PROGRAMS 

This section lists the messages printed by the programs and indi- 
cates what to do when they are encountered. The first group of messages 
(1 to 15) are concerned with errors in input quantities and are self- 
explanatory . 

(1) INTERVAL SIZE FOR PRESSURE DISTRIBUTION PRINT-OUT MUST BE GREATER 

THAN 0 AND LESS THAN 1 

This message indicates that the condition 0 < XLOUTP/i < 1 has 
been violated. 

(2) XMTB MUST BE GREATER THAN 0 AND LESS THAN 1 

This message indicates that the condition 0 < ( x /-^R max < 1 
has been violated. 

(3) EQUIVALENT BODY THICKNESS RATIO MUST BE GREATER THAN ZERO. 

This message indicates that the condition D / £ > 0 has been 
violated . 

(4) XMTW MUST BE GREATER THAN 0 AND LESS THAN 1 

This message indicates that the condition 0 < ^ ^ ) < 1 has 

been violated. w \ a x 

(5) WING THICKNESS/CHORD RATIO MUST BE GREATER THAN ZERO 

This message indicates that the condition Z /c > 0 has been 

max' w ^ 

violated . 

(6) XRLE MUST BE GREATER THAN 0 AND LESS THAN 1 

This message indicates that the condition 0 < X /£ < 1 has 
been violated. 

(7) CRT MUST BE GREATER THAN 0 AND LESS THAN 1 

This message indicates that the condition 0 < C /£ : < 1 has been 

K t 

violated . 
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(8) XRLE MUST BE LESS THAN XRTE 

This message indicates that the condition x r £ e < X rte has been 
violated. 

(9) ANGLE MUST BE BETWEEN 0 DEGREES AND 90 DEpREES 

This message indicates that the condition 0 < ANGLE < 90 has 
been violated. 

(10) AXIAL LOCATION OF BODY BASE MUST BE AT OR BEHIND POINT WHERE 

WING TRAILING EDGE PIERCES BODY SURFACE 

This message indicates that the condition X „ < X. has been 

rjZei d 

violated. 


(11) AXIAL LOCATION OF BODY BASE MUST BE AT OR BEHIND TRAILING EDGE 
OF WING TIP CHORD 


This message indicates that the condition 
violated. 


X sm 2 < *b 


has been 


(12) RATIO OF MAJOR TO MINOR AXIS MUST BE GREATER THAN 0 


This message indicates that the condition 
violated. 


X(=a/b) > 0 has been 


(13) TAPER RATIO MUST BE BETWEEN 0 AND 1 

This message indicates that the condition 0 £ TR £ 1 has been 
violated . 


The following error messages, numbers (14) through (18) should 
not occur in the present programs. If they do, they are probably caused 
by an error in reproducing the source decks. 

(14) EXECUTION TERMINATED BECAUSE EXPONENT N CANNOT BE DETERMINED 
TO WITHIN .01 PERCENT IN 20 ITERATIONS 


This message is printed by the circular body programs when the 
exponent n describing the equivalent body ordinates (see eqs . (11) to 

(14)) cannot be found accurately from information regarding the point of 
maximum radius (eqs. (12) or (16)) by using a Newton-Raphson 
procedure 20 times. 
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(15) EXECUTION TERMINATED BECAUSE POINT WHERE WING LEADING EDGE 

PIERCES BODY CANNOT BE FOUND TO WITHIN .01 PERCENT IN 20 ITERATIONS 

This message is printed when the point X . / g. cannot be found 

L 

accurately by using a Newton-Raphson iteration procedure 20 times. 

(16) EXECUTION TERMINATED BECAUSE POINT WHERE WING TRAILING EDGE 

PIERCES BODY CANNOT BE FOUND TO WITHIN .01 PERCENT IN 20 
ITERATIONS 

This message is printed when the point X . /I cannot be found 

rre 1 

accurately by using a Newton-Raphson iteration procedure 20 times. 

(17) EXECUTION TERMINATED BECAUSE SEB"(X)=0 POINT CANNOT BE DETERMINED 

TO WITHIN SUFFICIENT ACCURACY IN 10 ITERATIONS 

This message is printed by the circular body programs when the 
point x s /^ cannot be found accurately by using a Newton-Raphson iteration 
procedure 10 times. 

(18) INTEGRATION TERMINATED BECAUSE ACCUMULATED ERRORS HAVE CAUSED 

INTEGRATION SUBROUTINE TO BISECT ORIGINAL, STEP SIZE (.001) 10 
TIMES 

This message is printed when the integration subroutine used 
(Hammings modified predictor-corrector scheme as described in the 
Scientific Subroutine Package from IBM Corporation) cannot achieve the 
integration accuracy (DPRMT(4)) desired even though the original step 
size (DPRMT(3) = 0.001) has been bisected 10 times. 

(19) PROGRAM TERMINATED BECAUSE INDENTED BODY RADIUS HAS BECOME LESS 

THAN ZERO AT x/ g = 

This message is printed by the program when the cross-sectional 
area of the wing is larger than that of the equivalent body, so that 
the indented body radius is less than zero. A wing with smaller 
thickness/chord ratio or smaller span must be used. 

(20) PROGRAM TERMINATED BECAUSE INDENTED BODY MAJOR AXIS HAS BECOME 

LESS THAN ZERO AT x/f = 

This message is printed by the program for the same reason as 
message ( 20 ) . 
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(21) CALCULATION TERMINATED BECAUSE FLOW ABOUT EQUIVALENT BODY HAS 

BECOME SUPERSONIC AT x/ i = . INPUT MACH NUMBER GREATER THAN 

LOWER CRITICAL 

This message is printed when using the local linearization method 
for purely subsonic flow to calculate u B and indicates that the 
free stream Mach number is greater than the lower critical. There exists 
a region of supersonic flow and the local linearization method does not 
apply. 

(22) CALCULATION TERMINATED BECAUSE FLOW ABOUT EQUIVALENT BODY HAS 

BECOME SUBSONIC AT x/i = 

This message is printed using the local linearization method for 
purely supersonic flow to calculate u^ and it indicates that a 
region of subsonic flow has been encountered. This always occurs near 
the tail of the equivalent bodies considered here and the phenomenon is 
discussed fully in reference 6. However, if this region occurs at or 
near the nose, the local linearization method does not apply. In this 
case, the following error message is also printed, 

INPUT MACH NUMBER LESS THAN UPPER' CRITICAL 

indicating that the condition M < M has been violated. 

cr,u o° 

(23) START OF SUPERSONIC CALCULATION 

SUPERSONIC CALCULATION STARTS AT x/i = 

This message is printed when using the local linearization method 
for M^ » 1 flows to calculate u B and indicates where transfer is 
made from the parabolic (eq. (5) ) to the hyperbolic (eq. (7) ) differential 
equation. See reference 6 for details. 

(24) START OF SUBSONIC CALCULATION 

SUBSONIC CALCULATION STARTS AT x/i = 

This message is printed when using the local linearization method 
for M « 1 flows to calculate u D and indicates where transfer is 
made from the hyperbolic (eq. (7)) to the elliptic (eq. (6)) differential 
equation. See reference 6 for details. 
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NUMERICAL EXAMPLES 


General Description of the Output 

The output format of the program developed is as follows. On the 

top of the first page, a heading is printed describing the Mach number 

range, the general class of body (circular or elliptic) and wing being 

considered, and 'the theory used. Next, th^ wing-body geometry and 

flow field characteristics are printed. Then, if the local linearization 

method is being used, the program prints the fact that the integrations 

are starting at x = x g and proceeding to the nose. If the user 

supplies the distribution of u B , this is omitted. A heading of 

independent and dependent variables is printed next which contains, from 

left to right, the axial location x/£ at which output is to be given, 

the actual body radius R^/ £ (or actual semimajor axis a/ 2 . in the 

case of elliptic bodies) at that axial location, the angles 9 (in 

degrees) in the crossflow plane at which output is desired, the surface 

pressure coefficient C (body) , and six flow-field pressure coefficients 

P 

Cp(r/D = ) at the indicated distances r/D in the crossflow plane. 

For the case when the user supplies the u D distribution, the calculation 
begins at a point close to the nose and proceeds toward the body base 
with the values of the above quantities being printed out in the indi- 
cated tabular form at specified axial locations. When the local lineari- 
zation method is used to determine u D , the calculation begins at x = x 
and proceeds, to a point close to the nose, with the quantities described 
being printed at the specified axial locations, then, when the point 
close to the nose is reached, the program returns to x = x g , prints 
the fact that the integration is restarting at that point and proceeding' 
to the tail, prints the independent and dependent variable heading 
described above, and proceeds with the calculation to the body base. 

If it should happen at some point that the radial distance in the 
crossflow plane at which output is desired is less than the body radius 
(i.e. the point is inside the body), or if an output point falls on 
the wing leading or trailing edge, the pressure coefficient at that 
point is set equal to l.E + 6 and the program continues. Also, because 
of the discontinuity in the second derivative of the indented body 
d 2 R,/dx 2 (or d 2 a/dx 2 ) at the points x/£ = X . and x/f = X 
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and the discontinuity in the slope of the wing span at the points 
x/i = X sm , and x/i = X sm2 , output is not printed within a band 
(Ax/i = 0.005) of those points. When the calculations are. successfully 
completed, the pertinent aerodynamic coefficients are calculated and 
printed and the program then proceeds to read the data for the next 
case. 

Sample Cases 

In order to provide checks on the programs, sample test cases 
have been run for each program and the results are provided in figures 8 
through 12. In each case, the input data is provided together with 
the corresponding output. 

Finally, we note that in order to improve their accuracy, changes have 
been made in the subroutines, as given in reference 6, which compute the 
derivatives of the indented circular and elliptic body area distributions. 
Consequently, the test case results that appear in reference 6 will differ 
slightly from the results which the current program will produce for those 
same wing-body geometries. 
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Figure '2 i al . - General configuration of finite thickness 
wi eg- indented circular body combinations with 
st ra ight/swept forward trailing edge planforms 


Cross section for 



Figure 2(b) - General configuration of finite thickness 
wing-indented circular body combinations with 
sweptback trailing edge planforms 
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Figure 4.- Theoretical surface and flow field pressure distributions 
at = 1 for a noniifting parabolic-arc profile wing — indented 

parabolic-arc body combination; equivalent body thickness 
ratio D/f = 0.1, wing- aspect ratio AR = 1.7, thickness/ 
chord ratio t/c = 0.04, planform taper ratio TR = 0 . 2 , 

and with X r£e/j , = 0. 25, C = 0.50, X b/J 


0 . 86 . 






Figure 5.- Theoretical surface and flow field pr 
distributions and loadings at M = 1 and a 
for a parabolic-arc profile wing — indented 
parabolic-arc body combination; equivalent 
body thickness ratio D /i = 0.1, wing 
aspect ratio AR = 1.7, thickness/chord 
ratio t/ c w = 0.04, planform taper ratio 

TR = 0.2, and with * r fe/l = 0.25, 

d ~ 0.5, Xy. / 0 ~ 0.86. 








Figure 5.- Concluded. 
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LEGEND 



Wing-indented body, 9 - 
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Figure 6.- Theoretical surface and flow field pressure 
distributions $t = 1 for a nonlifting parabolic- 

arc profile wing--indented parabolic-arc body 
combination; equivalent body thickness ratio 
D /} - 0.1, wing aspect ratio AR = 2.8, 

thickness ratio t/c = 0.04, planform 

w 


taper ratio TR = 0.4, 


r 9 e / 


. = 0.25, 


and with 
0.3, 








Body surface 





LEGEND 

— — — — Wing-indented body, 9=0° 

— — Wing-indented body, 9 = 90° 

■ — — — Equivalent elliptic body, 9=0 
— Equivalent elliptic ' ody, 6=9 


Figure 7.- Theoretical surface and flow field pressure 
distributions at M = 1 for a nonlifting parabolic- 
arc profile wing — indented parabolic-arc body 
combination; having a body of elliptical 
cross section with ?> = 3 ; equivalent 
body thickness ratio D / f. = 0 1, wing 
aspect ratio AR = 2.8, thickness 
ratio t/c = 0.04, planform taper 

ratio TR = 0.4, and with 


'vie/ I 


= 0.25, C = 0.3, 


0 




-0. 



Figure 7.- Concluded. 
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fcTRAWIM AMAOJal. ,M0PT=1 , TAUB- . 1 ,TAU1^ . 04 , XHTB“ . 5 , XMTV** . 5 , AMGLE“58 . , 
SSMAX- . 25 , XRLE* . 2 5 , TR- . 2 f CRT- . 5 , XLBA6E- . 86 , XLOUTP- . 0 5 , 4 END 


(a) Input. 


CALCULAT1CN OF SURFACE AND FLOW FIELD PRESSURE DISTRIBUTIONS 
FUR FLOW AT FREE STREAM MACH NUMBERS AT OR NEAR ONE* BlLGW 
THE L CwER CRITICAL, OK ABOVE THE UPPER CRITICAL ABCUT A FINITE 
THICKNESS MING- INDENTEO CIRCULAR BODY COMBINATION hi ITH THE 
EQUIVALENT 000V OF REVOLUTION EITHER USER-SPEC IF |£0 OR HAVING 
URO IN AJ t S R PROPORTIONAL TO X/L-IX/L)**N OR 1-X /L- (1- X/L > **N . 
THt WING HAVING A CONSTANT T H I C XNE S S/CHCRO RATIC* TAPER RATIO 
BETWEEN C AND 1. ANC WITH 0R0 1 NATES l PROPORTIONAL TO XBAR/C 
-I X0AR/C l**N OR l-XBAR/C-il-XBAK/C l**M BY USING THE TRANSONIC 
EQUIVALENCE RUlE ANC THE LOCAL LINEARIZATION METHGO 


NING-BCOY COMBINATION CEOMETRY AND FLOW FIELD CHARACTERISTICS 


EQUIVALENT BUOY THICKNESS RATIO * 
EQUIVALENT BODY MAXIMUM THICKNESS AT X/L ■ 
EXPONENT N FOR EQUIVALENT BODY OROINATES » 
SEB* * ( X I • 0 U X/L - 
WING MAX. THICKNESS AT XBAR/C * 

WING THICKNESS/CHORD RATIO - 
EXPONENT M FUR WlNG URClNATbS - 
LEADING EDGE CF WING RCGT ChCRC AT X/L * 
TRAILING EOOE CF WING ROOT CHORE AT X/L - 
PLANFCRM TAPER RATIO » 

LEADING EDGE PlfPCES OCDY AT X/L * 

TRAILING A.OGE PIERCES eOCY AT X/L * 

BOOY RASE AT X/L - 

LEADING EDGE SWEEP ANGLE IOEGI - 

TRAILING EDGE SWEEP ANGLE (OEGl = 

LCCATICN OF WINGTIP LEADING EDGE AT X/L “ 
LOCATICN OF WINGTIP TRAILING EDGE AT X/L - 
NORMALIZED MAX. SEMISPAN SSMAX/L - 
ANGLE OF ATTACK ALPHA (OtG) - 
RATIC CF SPECIFIC HEATS » 

FREE STREAM MACH NUMBEP » 


O.IOOOOE 00 
0.50Q00E 00 
0. 200006 01 
C. 211326 00 
C.500COE 00 
C. *00006-01 
C. 20000c 01 
C.250C0E 00 
C. 1500C6 00 
C.20000E 00 
0.31960E 00 
C.75008E 00 
0.86000E 00 
C. 58QCGE 02 
O.COOOO 
C. 650006 00 
0.75008E 00 
C. 250006 00 
O.COOOO 
C.i*OOOE 01 
C.10000E 01 


START CF INTEGRATION FROM SE b • * I X ) - 0 TO NCSE 


X/L R BUOY /L THETAI0EG1 CPIBODY) CP(R/0« l.OOICPIR/O- 2.00ICPIR/D- 3.00ICPIR/0- A.OOICPIR/O- 5.00ICPIR/0- 6.001 


0.2113 

0.2113 

0.0333 

0.0333 

o.ccoo 
9.0C00E 01 

2 . 1 30 76-02 
2 • 1 307E- 02 

3. J159E-02 
3.J159E-02 

3. *2706-02 
3. *2706-02 

3 .**766-02 
3. **766-02 

3 .*5*86-02 
3.*5*8E-02 

3. *5026-02 
3. *5826-02 

3. *6006-02 
3. *6006-02 

0.2003 

0.2003 

C.C320 

0.0320 

c.occo 

9.0C00E 01 

2. Bbi 3E-02 
2. 88136-02 

3.81 716-02 
3.017 IE -02 

3.712*6-02 

3.712*6-02 

3 .60696-02 
3.6069E-02 

3.52*76-02 

3.52*76-02 

3. *5076-02 
3. *5076-02 

3.40396-02 

3.40396-02 

0.1503 

0.1503 

0.0255 
0.02 55 

C.0000 
9 • OCOJE 01 

6 . t692E-G2 
6. 66926-02 

5.9*716-02 

5.9*716-02 

*• 7*726-02 
*.7*72E-02 

* .007 16-02 
*.00716-02 

3. *7556-02 
3. *7556-02 

3. 06136-02 
3. 0613E-02 

2.72216-02 
2. 72216-02 

0.1C03 

0.1003 

0.01 81 
0.0181 

0.0000 
9 .OQOOl 01 

I.137*6-0i 

1.137*6-01 

7.5672E-02 
7.567 26 -02 

5.087*6-02 

5.087*6-02 

3.6118E-02 
3 .0 1106-02 

2.56C7E-02 
2.56C7E-0 2 

1.7**16-02 

1.7**16-02 

1.07656-02 

1.07656-02 

0.0*03 

0.0503 

U.C096 

0.00*6 

0.0000 
9.0C00E 01 

1. £1386-01 
1 • 31 3PE-01 

7 .9*6*2-02 
7.9*6*6-02 

*« 01 33E-02 
*.01 JJk-02 

1 .70366-02 
1.7038E-02 

6. 36126-0* 
6.3012E-O* 

-1 .2090E-02 
- 1 . 2090E-02 

-2.2*506-02 

-2.2*506-02 

6.00*3 

0.00*3 

0. C0C9 
0.0000 

o.ocoo 

9.0C00E 01 

3. 77*7E-CI 
3. 7 7*7t- 01 

* .6 1 C7E -02 
*.61076-02 

- 1.91176-03 
-7.9117E-03 

-3 .951 2E-02 
-3.9S12E-02 

-6 . 1932E-02 
-6.1932E-02 

-7. 93236-02 
-7.9323E-02 

-9.35336-02 

-9.35316-02 

START Lf 

INTEGRA T 1 

ICN FHCM SEO* 

•IX) ' C TC 

TAIL 






X/L 

P BODY /l 

fHL T A10EG) 

CPIbCOYI 

CPIR/D* 1.00)CP(K/0> 2.30ICPIR/0- 3.00ICPIR/D- *.00)CPIR/0» 5.001 

CP 1 R/0* 6.00) 

0.2113 

0.2113 

0.C311 
0.03 33 

0.0000 

9.0COOE 01 

2. 13076-02 
2.1307E-C2 

3.3159E-02 
3.31 59E-02 

3. *2706-02 
3. *2706-02 

3.**76E-02 
i. 4*766-02 

3. *5*06-02 
3. *5*86-02 

3. 4SE2E-C2 
3. *9826-02 

3.46006-02 

3.46006-02 

0.2503 

0.2503 

C.C3 15 
0.03^5 

c.ocoo 
9.0000E 01 

- 3. *5*5E-C 3 
-3.-5*5E-C3 

1 .A990E-C2 
1 .*9906-02 

2.30306-02 

2.3030E-02 

2 .731 16-02 
2 .731 IE-02 

3. 02786-02 
3.02 70E-O2 

3.25596-02 
3.25596-0 2 

3.4*136-02 
3. 4*136-02 

0.3003 

0.3003 

0.0*20 

0.0*20 

o.ocoo 
9.0C00E 01 

-3.165BE-C2 

-3.165BE-C2 

-8.32C7E-J3 
-8 ,3c 07E -03 

6.98*06-03 

6.98*08-03 

1 .5*006-0 2 
1 .56006-02 

2 . 16566-02 
2.1656E-02 

2 .63 37E-02 
2.63376-02 

3.015*6-02 

3.015*6-02 

0.3503 

0.3503 

3.0*53 
0. C*53 

0.0000 
9.00006 01 

- * . 2756 £- C2 
-9. 1*826-02 

-1.1*706-02 

-3.97366-02 

-*•*7726-03 

-1.131*6-02 

5.01*76-03 
1 .97686-03 

1.24596-02 

1.07*96-02 

1.04396-02 

1.73**6-02 

2.3*466-02 

2.26466-02 


Figure 8.- Sample input/output for a wing-body combination having a 
circular body and with TR f 0, B te 0, a - 0. 
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0*4003 

0.0463 

0.0000 


-1.3543E-C1 

-9.07596-03 

“2. 93876-02 

“l .64346-02 

-5. 040.6-03 

2.75471-03 

5.92106-0) 

0*4003 

0.0463 

9.00005 

01 

-1. 34216-01 

-7.92456-02 

-4.04 776-02 

-2.13096-02 

-8.575 TE-CJ 

1.00436-03 

8.70546-03 

0*4503 

0.0452 

o.ocoo 


-l. 65655-01 

-l .06616-01 

-4.7339E-02 

-3.4247E-02 

-2.22176-02 

-1.22 15E-02 

-3.93216-C3 

0.4503 

0.0*52 

4.00006 

01 

-1.6U5E-QI 

-I.Q676E-01 

-c. 3*796-02 

-4.LQME-02 

-2.60 16E-02 

- 1 .4? OlE-02 

-5.61536-03 

0.3003 

0.0421 

O.OCOO 


“2.02 766-0, 

-1.52426-01 

-5.A0S95-C2 

-4 .754ft -02 

-3. 58586-02 

-2.59C26-02 

-I. 74766-02 

0*3003 

0.0421 

R.OCOOc 

U1 

- 1.66336-01 

-1.18406-01 

-7.6063E-02 

-5 .65036 -02 

4.02 066-02 

-2 .86626-02 

-1.93R7E-02 


START OF SUPERSONIC CALCULATION 

SUPERSONIC CA iCUl A T I CN STARTS AT X/L - 0.542326 00 

X/L RBUOV/L THE T A t OfcG I CPIBCOVI CPIR/O* l.OOICPCR/D- 2.0QICPIR/D- 3.C01CPIR/0- 4.001CPTR/D- 5.00ICPIP/C- 6.001 


0.3503 

0.3503 

0.03 74 
0.C3 74 

0.0000 

9.0C00E 

01 

rl. 74066-01 
. -1.3590E-CI 

- 1 .60906-01 
-1.06336-01 

-4.**0*»7E-02 

-8.00846-02 

-5.5216E-02 

-6.19586-02 

-4.57295-0 2 
-4. 91 526-02 

-3. 72166-02 
-3. 93466-02 

-2.99535-02 

-3.14186-02 

0 .6003 
0.6003 

0.C319 

0.0319 

0.0000 
9 .00006 

01 

-7.e8466-02 

-4.49MIE-C2 

-1.33886-01 

-6.074OE-O2 

“1.03866-01 

-6.46226-02 

-5. 68455-02 
-5. £>7216-02 

-5.077 7E-C2 
-4.90856-02 

-4.49206-02 

-4.42336-02 

-I.95965-C2 

-3.94905-02 

0.7003 

0.7003 

O.OJOO 

0.0300 

O.OCOO 

9.00006 

01 

1.37146-01 

1.88036-01 

-6.96176-02 

2.08596-02 

-2.03766-01 
-3.81 87C-02 

-1.3049E-01 

-5.49145-02 

-9.94916-02 

-6.01426-02 

-8.61C3E-02 

-6.13455-02 

-1.80816-02 
-6. 10046-02 

0.8003 

0.8003 

0.0320 
0 • 0320 

0.0000 

9.0C006 

01 

81246-02 

-8.61246-02 

-7.39226-02 

-7.89226-02 

-8.0057E-82 
-8. 00576-02 

-8.122 7E -02 
-8.1227E-02 

-8.21 C2E-02 
-8.21C26-02 

-0. 28041-02 
-8.29C46-02 

-8 .33866-02 
-9. 3386E-02 

0.8503 

0.8303 

0. 0255 
0.0253 

O.OCOO 

9.0000E 

01 

-2.95356-C2 
-3. 95356-02 

-4 . 7040E-02 
-4.70^06-02 

-5.91916-02 

-5.919x6-02 

-o .6 66 OE -02 
-6.66806-02 

-1.20575-C2 

-7.2057E-02 

-1. 62475-02 
- 7. 62475-02 

-7.96796-02 
-T. 96 795-02 


WAG COEFFICIENT • 0.10440E 00 


(b) Output. 


Figure 8.- 


Concluded . 
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i TRAN IN AMACH=1. , *0PT«1 , TAL'B- . 1 , TAUW-. 04 , XMTB- . 5 , XMTW-. 5 , ANGLE-S0 . , 
SSKAX-. 25 , XRLEo. 2 $ , TR-. 2 , CRT- . 5 , X LEASE- .86, X LOU TP- . 05 , ALPHA-2 . , 6 END 


(a) Input. 


CALCULATE 01 SURFACE AN J FlllU FIELO PRESSURE DI $ 1 9 I BUT ION S 
FOR FLOW AT FREE STREAM MAC H NUMBfcH $ AT UP NEAR C.tot , BfcLCw 
THE LOWER CRITICAL, CK ABOVE THE UPPfcK C R 1 1 I CAL ABLLT A C INITE 
thickness v iso- indented circular buly ccmpinaticn with the 

ECU I V At E NT POUT OF REVOlUHCN EITHER USER - SPfc C I F I fc D OR HAVING 
JPOINATtS H PROPCkriCNAL TO X/L-IX/L)*** C« l -X/L - 1 l- X /L I **N * 
Tmc MING HAVING A CONSTANT THICKNESS/LHCRO RATIC, TAPER RATIO 
BETWEEN 0 AND 1, AND WITH ORDINATES Z PROPORTIONAL TO XBAR/C 
-lxeAR/C»**M CR l-XT)AR/C-li-X BAR/C i)Y USING T Ht TRANSONIC 

EQUIVALENCE RULE ANC Th£ LOCAL LINEARIZATION mETHCD 


WING-BODY COMBINATION OECMfclkV A NO HCh FlfctC CHARACTERISTICS 


EQUIVALENT BOOT THICKNESS RATIO * 
EQUIVALENT BODY MAXIMUM THICKNESS AT X/L • 
EXPONENT N F CR EQUIVALENT BCCV CHDlNATES «= 
5E8* M XI ■ 0 AT x/L ■ 

WING MAX. THICKNESS AT XHAk/C = 

WING THICKNESS/CHCRD RATIO * 

EXPONENT M FOR WING CMCINAUS * 

LEADING LOOT OF WING kCOl CHCRU AT x/L « 
TRAILING EDGE CF wING ROOT CHORD AT X/L « 
PLANFCRM TAPER RATIC * 

LEADING EDGE PIERCES RCOV AT X/L - 
TRAILING EDGE PIERCES eCCV AT X / 1 « 

BODY BASE AT X/L » 

LEADING EDGE SwECP ANGLE ICfcGl * 

TRAILING EDGE SwEEP ANGLE (OtGl - 
LCCAT1CN OF WINGTIP L fc AC I NG EDGE AT X/L - 
LOCATION OF WINGTIP TRAILING ECCL AT X/L * 
NORMALIZED Max. SF.M1$PAN SSMAX/L » 

ANGLE CF ATTACK alpha (DEG I - 
RATIC CF SPECIFIC HEATS ■ 

FREE STREAM MACH NUMBER * 


C.lOOOCfc 00 
0.500GQE 00 
C.20000E 01 
C. 211326 00 
C . 500 OUE 00 
C.40000E-01 
C . 200QCfc 02 
c. 250006 00 
C.TSGUCE 00 
0.2000CE OC 
C. 319606 00 
C. 750086 00 
C.8N00CE 00 
C.5B0CCE 02 
C. 00000 
C.6500BC 00 
C.75009E 00 
C • 250C0E 00 
C.20000E 01 
0.140GQE 01 
0.10000E 01 


START OF INUGRATICN MUCH SEP'MXI « C TC NCSE 


X/L 

RBOOV/L 

TMETAlOfcGl 

CPTSCCYI 

CPIR/O* | .00 I C P | R/()« 2.0U)CPIk/U« 3.C0ICPIR/D- 4.0CICPCR/0* 5.00ICPIM/0- 6.00 

0.2113 

O.C333 

o.ocoo 


1. 7652 6-02 

3 . 28766-02 

3. 42026*02 

3.4446E-02 

3.4531E-C2 

3. 457IE-C2 

3.4502E-U2 

0.2113 

0.0333 

5 . OCOOE 

01 

6. 40316-03 

2 . 56 5 2c -02 

3.U344E-02 

3.1830E-02 

3. 25546-02 

3.29B2E-02 

3 . 32t5fi-02 

0.211 i 

0.0333 

-** . OCOOE 

01 

2 • 8669 E-C2 

4.11 78E-02 

3.8321E-C2 

3.71826-0? 

3.6576E-C2 

3.62C2E-02 

3.59496-02 

0.2003 

0.O320 

o.ocuu 


2.5157E-C2 

3 . 79C8E- 02 

3.7O0OE-O2 

3 .o 04 lfc-02 

3.5231E-02 

3.4577E-02 

3.40 J2E-02 

0.2003 

0.0320 

9.0000E 

01 

1.3296E-C2 

3.0660E-02 

3. 31 98C-02 

3.3426E-02 

3.3256E-02 

3.299JE-02 

3.27006-02 

0.2003 

0.0320 

-9.0CCOE 

01 

6.67606*02 

4.61 77E-02 

4.11736*02 

3. 87686*02 

3.7269E-02 

3.62C4E-02 

3.53856-02 

0.150 » 

C. 0255 

C .OCOO 


6. 30>7E-C2 

5.S3C7E-02 

4. 74 326*02 

4 .005 i£-0 2 

3.4745fc-02 

3.06C7E-C2 

2.72166-02 

0.1503 

O.U255 

0 . OCO )t 

01 

4.8381E-02 

5.2306E-02 

4.37S0E-O2 

3 .7 600E-02 

j . 2B97E-02 

2.91246-02 

2.59796-02 

0.1503 

0.C255 

-9. OCOOE 

01 

8. 764 IF.- C2 

6 .6 9456 -02 

5. 12 33 E* 02 

4.2577E-02 

3.6634E-02 

J.2USE-02 

2.9472F-02 

0.1003 

0.01 81 

U.OUOO 


1 • i OOB fc- 01 

7. 55*16-02 

5. 0<)54t-02 

J .o 1 Q9£-02 

2.56C2E-02 

1.74386-02 

1.0762F-02 

0.1003 

o.ci ei 

4.0CCOE 

01 

9.26*1 E-02 

6 .97 7 IE -02 

4. 7879E-02 

3.4 11 *E-02 

2.41 C2E-02 

1.6236F-0? 

9. 75S7E-0 3 

0.1003 

0.01 B1 

-9.UOOOE 

01 

1.372RE-C1 

8. 172M L- - 02 

5. J9C7c- 02 

3 • 39E-02 

2.712^6*02 

1 . 8653 k - 02 

1<1 774E-02 

0.U5C1 

0. C096 

0.0000 


1.77736*01 

7.9642C-02 

4.01 2bc*02 

l. 70356-02 

6. 3473E-C4 

-l. 20916-02 

-2.2451E-02 

0.0503 

0. 0096 

9.0000F 

01 

1.5/686-01 

7.5897t-02 

3.B340E-02 

1.584CE-02 

-/.6/43E-04 

-1.2809F-02 

-2.3090E-02 

0.0503 

0.C0S6 

-9. OCOOE 

01 

2. 0771 fc* Cl 

0.3O7fcE-O2 

4. 193PE-02 

1.F240E-02 

1.5375E-03 

- 1. 13696-02 

-2. 10SCE-U2 

0 .006 3 

O.COC9 

U.OCOO 


2 . 738 l£*Cl 

4 ,61 C7F-0 2 

* 7. 9l 1 7E-0 J 

*3.9512 c -02 

-6.19326-02 

- 7.9323E-02 

-9.35335-C2 

0.0C-* 3 

0.0009 

9. OCOOE 

01 

3.5IU0E-01 

4 .57506 -02 

*8 . 09 C2 E-Oi 

*3.963 16*02 

-6.2022E-02 

-7. 93956-02 

-9. 35S2 r -02 

0.0063 

0. C0C9 

-9.0000E 

01 

6.063 7fc- C 1 

4 .64636 “02 

-7. 7331 6-03 

-3. 93936*02 

-6. 1843E-02 

-7.9252E-C2 

-9.3473C-02 

START OF 

INTEGRATION fmcm 

SL B * 

•1X1 « C TC 

TAIL 






X/L 

RitlDY/L 

1 »(f T A IP 

tG 1 

OPIl.l CYI 

CP (R/n« 1 * Out f P |R/C» 2 . 00 ) L P ( k / D * 3.00)CPIR/D> 4.O0ICPI R/0- 5.00KPIR/D* 6.00 

0.2113 

0.C33J 

c.ccco 


l . 7 6 2 F - 0 2 

3. 28746-02 

3 . 4202 £-02 

3.44466-02 

J.4511E-02 

3. 45716-0* 

3.45S2F-C2 

0.2113 

0.0333 

9.CC00E 

O’ 

*• . 603 1 E-0.* 

2 .565 2t-02 

3.03 **4c*C2 

i.i 8*06-02 

3. 25546-02 

3 . 29626 - C2 

3 .32651 -02 

0.211 3 

0. 0333 

-9.0C0JI 

01 

2.6665E-C2 

4. H/8C -02 

3.833IL-02 

3.7 IB2E-02 

3. 65766-02 

3 . 6202F-02 

3.5949F-02 

0.2503 

0.0*75 

i . OCOO 


- /.liH9f-C3 

i . *• 6 2 3 E - >i 2 

< .2»42t-C2 

2.72 736*0/ 

3.0257F-02 

3.25456*02 

3.4w04f-02 

0*2503 

O.C I 75 

9 .0000 E 

0 1 

*l.61alE*C2 

7. 82/6E-01 

1.9/ *5 6-0/. 

2. 4 7466-02 

2.83436-02 

3.10056-02 

3.3116E-Q2 

0.2503 

O.C37S 

-9. OCOOE 

.1 i 

1.17udt-C2 

2.2? »lc-02 

/ . &9 S3 fc-02 

2.99ii£-02 

3.2257E-02 

3.* 1 40E-C2 

1 . 5 7306-02 


Figure 9.- Sample input/output for a wing-body combination having a 
circular body and with TR f 0 , B^ e < 0, a ? 0. 
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0.1003 

0.30<U 

0.1003 

0.U420 

0.0420 

0.0420 

o.ocoo 
9, OCOQE 01 

-9.ocooe oi 

-3.S3146-C2 
*4. 15426-02 
- 1 .S288C-C2 

-8.70916-03 
-1.45456-02 
-l. 11216-03 

6.97486-03 

3.62676-0) 

1.05936-02 

1 .5 55 26-02 
1.33196-02 
1.79756-02 

2.16296-02 
l . 9932(6—02 
2.34346-02 

2.63206-02 

2.49516-02 

2.77576-02 

3.01426-02 

2.09976-02 

1.13)66-02 

0.3509 

0.3503 

0.3503 

0.3503 

0.C453 
0.0453 
C.C4 53 
0.0*53 

0.0IUPP6R) 
0.0ILD4ER) 
9.CCOOE 01 
-9.0000c 01 

-1.C355C-C1 

2.04 72t-C2 
-1.15506-01 
-t.SQ25t.-C2 

-1.23156-02 
-1.23156-02 
-5 .62216-02 
-i. 02656-02 

-4.4)576-04 
- 4. t»3 976-04 
-2.21685-02 
-1.71506-04 

4.94516-03 

4.94516-03 

-5.52046-03 

9.60576-03 

1.24206- C2 

1.24206- 02 
5. C5C66-03 
1.65236-02 

1.84146-02 

1.84146-02 

1.27556-C2 

2.1961E-C2 

2.33896-02 
2.33896-02 
1. 0&C9E-O2 
2.65186 -02 

C.4QC3 
0.4003 
0.4003 
0.400 J 

0.C463 
0.0469 
0.0463 
0.04 63 

C.OltPPEK) 

U.OtLOMEA) 

q.ocoofc oi 

-9.00006 01 

-1.60536-01 

-7.99)06-02 

-1.65946-01 

-i.0004£-Ci 

-1.78006-02 
-1 .7800E-02 
-! .06476-01 
-5.07036-02 

-2.974*6-02 
-2.97446-02 
-5. 80796-02 
-2.21996-02 

-1.65776-02 

-l.*5776-02 

-3.38116-02 

-8.57426-03 

-5.915C6-C3 
-5.91 506-03 
-1.81946-C2 
1 .16816-03 

2.70776-C) 
2.70716-0) 
-6.7745E-C) 
8. 87146-03 

9. 88P66-03 
9.88866-03 
2.17926-0) 
1.52536-02 

0.4503 

0.4501 

0.4303 

0.4503 

0.C452 

0.0452 

0.0452 

0.0452 

c.ouppm 
O.OUUUER) 
9.0C006 01 

-t.ooeoE oi 

-2.3oC6t-Cl 
-l.))i8E-0l 
-1.95736-01 
- 1*241 46-01 

-1.80346-01 
-3 .39056-02 
-1.38256-01 
-7.3690E-02 

-4.812 7 t-OZ 
-4* Bl 276-02 
-8.57766-02 
-4.06756-02 

-4.45076-02 

-3.45076-02 

-5.73546-02 

-2.44006-02 

-2.235I6-C2 

-2.23516-02 

-3.87176-02 

-1.30946-02 

-l.2356t-C2 

-1.2)586-02 

-2.50586-02 

-4.19776-03 

- 9, 46846-03 
-1.48646-03 
-1.43416-02 
3.21376-03 

0.90CJ 

0.5003 

0.9003 

0.500) 

0.0421 

0.0421 

0.0421 

0.0421 

O.OItPPgR) 
0.0110*6*1 
9.0 COOC 01 
-9.0C00E 01 

-2.51)16-0 
-1.51 776-01 
-2. C2506-CI 
- l.27t>3f-0l 

-2.11306-01 
-9.21986-02 
-1 .52456-01 
-8.25736-02 

-5.79766-02 

-5.79766-02 

-1.0)756-01 

-5.144)6-02 

-4 .79906-0 2 
-4.79906-02 
-7 .49916-02 
-1.54906-02 

-3.60796-02 

-3.60796-02 

-5.56416-02 

-2.44466-02 

-2.40356-02 

-2.60356-02 

-4.13566-02 

-1.57506-02 

-1.75456-02 

-1.75656-02 

-3.01356-02 

-0.48376-03 


*f4#f Cf StPfft SO* K CAUUU TIC* 

$UP£MSCN1C CALC'Jl A T I CN STARTS AT X/L • 0.5*2326 00 


X/L RBJDY/L ThfTA(OEG) CPIBCCY) CPIR/D- 1.00>CPfft/Q» 2.Q0ICPIR/D* 3.00ICP ««/D» 4.00>CPf ft/0* 5,00>CPU*/D» 6.00) 


0 .5533 

0.0)74 

O.OlOPOCKt 

-2.2U66-C1 

-2.14716-01 

-5.3007fcr02 

-5.i*499fe-02 

-4.60736-02 

-1.74166-02 

-3.00856-02 

0.550) 

0.0)74 

C.OtLCtecRI 

-1.24536-01 

-1.05596-01 

-5.30C7 t-02 

-5.59946-02 

-<*.60736-02 

-3.74166-02 

-3.0085E-02 

0.5503 

0.0374 

9.00006 01 

-1.7)616-01 

-1.44186-01 

-l. 08516-01 

-8.42076-02 

-6.70676-02 

-5.42176-02 

-4.40796-02 

0.5503 

0.0)74 

-9.0C00E 01 

-4 • 57*36-02 

-7.05516-02 

-5.0508E-02 

-3.90206-02 

-3.07966-C2 

-2.41 74E-C2 

-1.8540E-02 

0.600) 

O.C319 

0.0IUMP6P » 

-l.22)5c-01 

-1.84246- 01 

-2. 22196-01 

-5.82336-02 

-5.12986-02 

-4.52046-02 

-4.01836-02 

0.6003 

0.014 

0.0<LGb£ft t 

-2.89046-02 

-8.16856-02 

6. 07 15 €-04 

-5.02336-02 

-5.12986-02 

-4.5ZC46-02 

-4.01836-02 

0.6003 

0.0)19 

9.CC00E 01 

-0,42876-02 

-1.06C7E-01 

-9.5)026-02 

-8.14006-02 

-7,00016-02 

-6. 11586-02 

-5.39866-02 

0.600) 

0.0)19 

-9.0C0Q6 01 

- 2 . 2376 6-C3 

-2.9)m-02 

-3.2bl2€-02 

-3.1 L95L-02 

-2.51276-02 

-2.44156-C2 

-2.4TC76-02 

0.700) 

0.0)00 

O.C(UPP£«7 

1.4021 6-Ci 

-6.7 8906-02 

-2.07)96-01 

-1.33266-01 

-1.00276-01 

-8.65096-02 

-7.8 ) 376 - 02 

0.7003 

0.0300 

o.o<n«e«» 

1.36516-01 

-6 *.935 76-02 

-2.0202E-01 

-1.33266-01 

-1.00276-01 

- 8 .45056-02 

-7.8)376-02 

0.7003 

0.0300 

9.0COU6 01 

1.67406-01 

2.16386-02 

-3.71936-02 

-5.40656-02 

-5.94486-02 

-6.0TS96-02 

-6.05376-02 

0.7003 

C.C3C0 

-9.0C006 01 

1.4110E-01 

2. 21866-02 

-3.7*436-02 

-5.47746-02 

-6. 01516-02 

-6.14456-02 

-6. 11146-02 


MAC CC66PIC16A7 - 0.134196 CO 
tiff CCEfnClBh T - 0.170706 01 
PITCH!*; HOHCHT CUffflCleWT * *0.090566 CO 


(b) Output. 


J 
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Figure 9.- Concluded. 



t TRAN IN AMACH° 1 . , MOPT-1 , TAUHa . 1 ,TAUW=. 04 , XMTB= . 5 , XMTW- . 5 , ANCLE-58 , , 

SSMAX- . 25 . XRLE- .25, TR- . 2 , CRT» . 5 , XLBASE- . 86 , XLOUTP- .05 .ALPHA-2 . , AL-3 . , (END 


(a) Input. 


CALCIUM ICr. ill- SURFAC r : AND FLOW FltLO PRcSSUKt ClilRlBUlIUNS 
FCR FLC* AT FKt£ S T R c AM MACH NU*B£»S AT OP SEAR DNE , BELOW 
THE LCwFP CRITICAL, UR ABOVt THE UPPER CRITICAL ArtCUT *. FINITE 
THICKNESS m INC- INCcN TEC RUDY f LjMb I SA T IUS WITH THE BCCY HAVING 
A‘: ELLIPTIC CROSS SECTION THAI MAINTAINS A CUNSIANT RATIC OF 
MAJCR/ M INCR AXES AL rjfiG TM C ^ N T I Rt BC.)Y LENGTH WITH ThE f.QUl- 

VALENI BCCV OF revoluticn cither user-specified cr havisc CR- 
O/NlTtS R PROPORTICNAL TO X/L-1X/L I *•« 0« l-X/t- < i-X/L )**N, 

THc wISC HAVING A CONSTANT TM ICkNESS/CHC»G RATIC, TAPER PATIO 
HeT rnteu 0 AN J 2 , A MU mllH ORDINATES l MQPCRT HJNAL TO X9AR/C 
-lXEAK/cl**K CR 1-XBAR/C-li - XBAR/C I • • H H Y USING THE TRANSONIC 
EQUIVALENCE RULE AND THE LOCAL LINEARIZATION **E THCD 


HI NG-bCDY COMBINATION ('.fc 0*6 TRY ANU FLOW FltlC CHAR AC TE R | S T 1C S 


RATIC CF SCMIMAJOw/SEM (K 1NCR AXIS * 
EQUIVALENT BODY THICKNESS RAT |C » 
EQUIVALENT eocv ‘'AXIHU"’ ThICKNFSS at x/l « 
EXPONENT N FOR EQUIVALENT BUCY CRUINATES « 
SC 8* * T XT * o AT x/t * 
wing rax. thickness at xbar/c « 

WING THICKNC SS/CHCRO RATIO » 

EXPONENT M FOR WING ONCINATES « 

LEADING F0C£ CF WlNf. "CUT CM2RC AT X/L * 
TRAILING EDGE CF MING ROOT CHCRC AT X/L » 

planfckk Taper ratio * 

LEADING cOCE PIERCES RCCV AT X/L * 

TRAILING I JOE. PIERCES EOCY AT X/L - 
SCO Y BASE AT X/L » 

LEAOING code SWEEP ANGLE I CEO T » 

TRAILING ecct SWEEP angle (OtGI • 

LCCATICN CF m I NG T 1 P LEADING ECCc AT X/l * 
LOCATICN OP WINGTIP TRAILING FCGfc AT X/L * 
NURMAL I ZLO max. se **i span SSMAX/L - 
ANGLE OF ATTACK ALPHA IJEG1 » 

RATIC CF SHCIFIC FEAIS - 
FREE STREAM MACH NC m 8»:R 5 


C.30000E 

01 

0. lOOOCE 

00 

C.SOOOCE 

00 

C. 20000c 

01 

C.2H32E 

00 

C. 5000CE 

00 

C.'.UOQOfc- 

-01 

C.2U0C06 

01 

C.25CCCE 

00 

C.75Q0CE 

00 

C.2J000E 

00 

C.3B07CE 

00 

C. 75 00 BE 

00 

c.esoccE 

00 

C.5H0QCE 

02 

C.COOOO 


C.65008C 

00 

C. 7500AE 

00 

C.250CCE 

00 

C. 20000c 

01 

C.lACOOfc 

01 

C.10000E 

01 


START CF INTtGRATlCN FRi:* Sfc9"<XI = 0 TC NCSE 


X/L RrtOCV/L THFTAIDEG1 


CPIBCJVI CPIR/D=* UOOICPIR/O* 2.00TCPIR/U" 3.00TCPIR/D- 9.00ICPIR/0* 5.00ICPIR/D- 6.001 


0.2113 

0.0577 

O.UOQO 


2.9c 4 7 6- 02 

3 .66546-02 

3.5139E-02 

3.9M58E-02 

3.9762E-02 

3.9718E-02 

3.9695E-02 

0.2113 

0.C1G2 

9.0CCUE 

01 

-9.99015-03 

l .92676-02 

2.7C71E-02 

2.97266-02 

3 • 1020E-C2 

3. 17B0C-C2 

1.22706-02 

0.2113 

C.0192 

-s.ocooe 

01 

J . 22 59t- C2 

9 . 202 7c -02 

3.98 G9E-02 

3.C996E-02 

3.76556-02' 

3.7119E-02 

3.67 336-02 

0.2QGJ 

0.C555 

C.OCOO 


3.EU27E-C2 

5 .2 12 16 -02 

3. B0S3E-02 

3.69786-02 

3.59 766-C2 

3.97336-C2 

3. 91906-02 

0.2003 

O.OIB5 

9.0C0OE 

01 

6. 74 JOE— G— 

2 • 392oE -02 

2. 9 0 55c- 02 

3 .1299c -02 

3. 17CSE-02 

3 • 1781F-02 

3.17156-02 

0.2C03 

C. Cl c5 

-9 . OCOOLi 

01 

3.53326-02 

9.6 07 3E -02 

9.26986-02 

9.0C6UE-02 

3* 8336E- 02 

3.7 IC56-02 

3.61626-02 

0.1503 

J.C992 

C.OCOO 


B. *,B 5f E-C2 

6.36C-.E-02 

9. ew 6 1 E-02 

‘,.050 76-02 

3. 5000E-02 

3.0765E-02 

2 .73296-02 

0.1503 

0.C197 

9.0C006 

01 

’.So? o:-C2 

9 . -9v)Ot-02 

s.09 716-02 

3.5539E-02 

3. 19 1BE-02 

2.7975E-C2 

2.5091E-02 

0.1503 

0.01-7 

-9.0CC0E 

01 

». 97116-02 

6. 73dSc-Oi 

5.26C3E-02 

9 .2 75 1£ -02 

3. 76 10E-02 

3.29916-02 

2.9185C-02 

0.1003 

0. 03 l 3 

C.OCOO 


l.Hl 30E-01 

7 ,d6|7t -02 

5. l->93E-02 

3 .69 3 7t -02 

2.57866-02 

l • 7556F-02 

1 .0899F-Q2 

0.1003 

0.010- 

9.OC00F 

Oi 

d.fc8».3t-02 

6. 396 ic-02 

9. 52C6l-02 

3.29066 -02 

A. 25206-02 

1.53176-02 

9.00916-03 

0.10U3 

0 . cic- 

-9.9CC01* 

01 

: . :b-ie-ci 

6 • 2So-»L- -02 

5.51 90E-02 

3.9 190F-02 

2.7535F.-C2 

1 • 5338E-C2 

1 .23626-02 

0.0503 

0. C 1 66 

O.OCJU 


2. 1 96 1 fc- C i 

B.C5O5E-02 

^ • C*, C7L- C2 

1.7 155E-02 

7.09336-09 

-1 .2097E-C2 

-2 .29606-02 

0.0503 

-i. 0055 

9.UOOOE 

3 1 

I.21B2E-C1 

7.2-.84E-02 

3.tH7A£-02 

l .95206-02 

-5.30 70E-09 

-1. 3333E-02 

-2.3570F-C2 

0.0503 

0.00*5 

-9.OCC0E 

01 

I . d* o 2t-01 

B.9>0«,E-02 

9./651F-02 

1 . 89 1 36 -0 c 

i . 0667E -Ci 

-1.09356-02 

-2.15216-02 

0 .OCA 3 

0. CO 15 

C.OCCC 


w.2o6Jfc-01 

9.61196-02 

-7. i0b4t-0? 

-3.-5106-02 

-6.15326-02 

-7 • 91 23E-02 

- 9 • 35 326-C2 

0 .005 3 

O.OOC5 

4.00001 

Oi 

J ■ C9^ Vf-Ol 

9.55016-0^ 

-6.2121L-03 

-3. c 7i iE-02 

-6.20926-02 

-7. 99936-C2 

-9.3632E-C? 

0.00-3 

J.0005 

-V.OCOOf 

.01 

3. 73.3 1?-C1 

9.66916-32 

-7.61 706-03 

-3.^3196-02 

-d. 17896-02 

- 7.92C9E-02 

-9. 39 39F-C2 


i I AR T CF 
X/L 

INTtCRAT 

*dn::v/i. 

I CN FFCW iC5«, 
T Ht TAOfcG | 

•4X1 * C TC 
CPldtCVI 

TAIL 

CPIR/D- 1 .U0TCP1R/L- 2.001CPIK/0- 3.CCTCPIR/0- 9.00»rP(S/C- S.OOICPIA/D- 6.00 

0.211 3 

C. U6 77 

C.OCOO 


2. 56 *76-02 

1.68546-02 

j.Sl 496-02 

3 .-.BSSE-C? 

3.97626-02 

1.97 1 6E-C2 

3.96556-02 

0.2113 

J. 0 1 t2 

G.GCCOt 

01 

-9.9oGlt-C * 

1 . 42d 76-02 

2.70 fir-02 

2.9726fc-02 

3.10206-02 

3.1780E-C2 

3.22786-0^ 

0.2113 

C.01 52 

-G.OCCJt 

01 

"2.2259F-C2 

9 . 202 76-02 

3.9B456-C2 

3.P596E-02 

3.765S6-C2 

3.7U9E-C2 

5.673JE-02 

0.250 3 

0.0650 

Q . C COO 


- 3. 338 5E- C9 

1.7161c -02 

2.35696-02 

2.75976-02 

3.09 116-02 

3.26996-02 

3.9972F-C2 

0.2503 

0.0217 

9. OCOOC 

01 

-2.J538C-02 

>.05616-03 

1 .63556-02 

2 .2 893E-02 

2.69256-02 

2.98 896-C2 

1.2190E-02 

0.2503 

0.021 7 

-5.CCC3? 

01 

d • 66546 - C3 

2.991 16-02 

2.87916-02 

3.13376-02 

3. 33 736-C2 

3.50681-02 

1.6522E -02 


Figure 10.- Sample input/output for a wing-body combination having an 
elliptical body and with TR + 0, B t# < 0, a * 0. 
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0.30C3 

C.C720 

O.OCCO 

-3.4l76t~C2 

-9. 35416-03 

4. 88 10 £-03 

1 .55616-02 

2.16356-02 

2.43246-C2 

3.01 *5E-02 

0,3003 

0.0243 

9 . OOOQt OJ 

-4. *5696-02 

-1 .66846- 02 

1.6510E-03 

l. 186 56-02 

l • 80C46-G2 

2.40386-02 

2.82276-02 

0.3003 

C.C24J 

-9. OCOOE Oi 

-1.0815f-C2 

l .6 1 996-0 3 

t • 25536-02 

1.94226-02 

2. 43536-02 

2.66466-02 

3.2I0OF-O2 

0.3503 

0. C7ftfi 

Q .0000 

-6 . J4A3E-02 

-3.68k5c-02 

-i,. 06851-02 

1.78346-03 

1.03726-02 

l. 69866-02 

2.2366E-02 

0.3503 

0.C263 

S. OCOOE 01 

-»».3144£-C2 

-3.499S6-02 

- 1 . 34656-02 

-5.9450E-O4 

8.4323E-03 

1 • 53576-02 

2.0967E-02 

0.3 503 

0.0263 

-9. OCOOE 01 

-*.3 8446-12 

-2.06J36-02 

-4, 71576-Oi 

5.50716-03 

1.3085E-C2 

1.9IC96-02 

2.41 C76-02 

0.4 00 3 

0 . 08 30 

0. J (UPPtK J 

- 1 .33 J8C-CI 

-7 , 78496-03 

-1.99*96-02 

-8.TQ76t-03 

2 • 37046-04 

7.41496-03 

1.33 7GC-C2 

0.4003 

0.0830 

O.OUOUEK » 

4.41 716-02 

-7.7*496-01 

-1.99o9E-02 

-8.70766-03 

2.37046-04 

7.4145fc-03 

1.3370C-O2 

0.4C0J 

O.C2T7 

S.CCGCE 01 

-l. 35186-C1 

-7 .8S&3E-02 

-4. 1570E-02 

-2.21106-02 

-5.4 JQ06-03 

-l. 16666-04 

7.21 24E-03 

0.4003 

0.02 77 

-9.0000E 01 

-7. 0O3CS-O2 

-3.28501-02 

-1.22686-02 

-J .35396-03 

6.50476-03 

1.27716-02 

1. 80186-02 

0.4503 

0.0034 

O.OtilPPER 1 

-2.24d4f-Cl 

- 1 «9656c - 0 l 

-5. J03PE-02 

-J.755U-02 

-2.44926-02 

-i . 34216-C2 

-5.12656-03 

'0.4503 

0.0836 

0 .Of LOWER ) 

-9.54296-02 

-4.741 3E-02 

-5. 3038E-02 

-3 .75516-02 

-2.44926-02 

-1.39216-02 

-4.12656-03 

0.4503 

C.C275 

s.cccae oi 

-1 .S837E-C1 

-l .35476-01 

-8.O272E-02 

-5.82456-02 

-3.95446-02 

-2.? 7156-C2 

-1.40C4E-C2 

0.4503 

0.02 >9 

-9, OCOOE 01 

-1.25386-01 

“7 .25526-02 

-4.2515E-02 

-2.62 7 IE -02 

-1.46806-02 

-5.47C76-03 

2. 23176-03 

0.500 3 

o.oaot 

Q.OIUPP IR ) 

-2.e>T87£-Cl 

-2. 39756-01 

-7.29126-02 

-5.89206-02 

-4.43166-02 

- 3. 22416-02 

-2.21336-02 

0.5003 

0.06CI 

a.oaouERi 

-1.50226-01 

-l.U36fc-Oi 

-7. 2912E-02 

-5.89206-02 

-4.43166-02 

-3.22416-02 

-2. 213 JE -02 

0.5C03 

0.C267 

5. OCOOE 01 

- 2. 3275E-CI 

-V .6946d“0l 

-1. U77E-01 

-0.49606-02 

-6 . 32806-C2 

- *i . 71*26-02 

- 3. 4 3S5F-C2 

0*5003 

0.0267 

-9. OCodE 01 

-1 . 55o36-Ol 

-1.00226-01 

-6.5ioee-o2 

-4.S9356-02 

-3.24696-02 

-2.1853E-C2 

- i . 30C3 r :-02 


SfART Cf SUP6RS0MC CAlCClAT tCN 

SUP6RSCNIC CAIXUIATICN STARTS At X/l * 0.54232E 00 


X/l 

RQQDY/l 

7H6 7 A 1 D6G > 

CP18CCY1 

0.5503 

0.0729 

0.0«JPPER 1 

-2.44836-01 

0.5503 

C.C729 

0.0 (LOWER! 

-1. J4J6t-Cl 

0.5503 

0.0243 

9 .00006 01 

-2. 18796-C1 

0.5503 

0.0243 

-9. OCOO£ 01 

-i • 39 J8C- Oi 

0.6003 

0.0631 

0.0 (UPPER I 

-1.3192E-C1 

0.6003 

O.ObJl 

O.Oa-lwER 1 

-2.019 IE-02 

0.6003 

0.0210 

9 • OCOOt 01 

-1.19926-Ci 

0.600J 

C. 0210 

-9.0C0C6 Oi 

-2.0lo3£-C2 

0.7 CO 1 

0.C550 

0.0(CPPE« > 

1.323QE-CI 

0.7003 

0.0550 

0.0 (low ER » 

l .268 96-Cl 

0.7003 

0.01 83 

9. OCOOE Ql 

1.5910E-C1 

0.7C03 

0.0183 

-9.0CCOE 01 

1.607U-01 

CRAG COEFFICIENT 

m 

C. 124626 CO 

1 1 FT CCFFFIC1ENT 

» 

C. 170736 Cl 

PITCHING 

PCMLNT 

CLcFFJClENT * 

-C.8f94CE CO 


CPIR/O* 1.00)CPtR/C- 2.00llP(R/0» J.CO 


-2.3So3c-OL 

-6. 74486-02 

-6.7246E-02 

- 1 .21**26-Q1 

-t. 74986-02 

-6.724t>E-02 

-1 .6 8786-01 

-1.2'»70E-0i 

-9.60746-02 

-9. 56&*E-02 

-6. 7l59t-C2 

-S.1242E-02 

- 1 . 7966E-0 l 

-2. 23966-01 

-5.915TE-Q2 

-7.52536-02 

-2.47&2E-04 

-5.51576-02 

-1.16826-01 

-5.97846-02 

-8,413 IE-02 

-4. 06306-02 

-3.7462E-02 

-3.42076-02 


-J.90G26-Q2 

-1.S151 6-01 

-1.24896-01 

-4.09516-02 

-1.86716-01 

-1.2409E-O1 

2-. 15666-02 

- 3» 1702 6-02 

-4.89176-02 

2.18036-02 

-3.22 10E-02 

-4.96256-02 


>fP|R/0* 4.0GICP<R/C- 5.00»CP(»/C« fc.OOt 


-5,47186-02 

-4.39B36-Q2 

-3.49326-02 

-5.47 lflfc-02 

-4,39836-02 

-3.49326-02 

-7.601 96-02 

-6.09566-02 

-4.90276-02 

-4.00276-02 

-3.11436-02 

-2. 3ft 826-02 

-5.22C56-02 

-4.5961E-02 

-4.07526-02 

-5.22C56-02 

-4.5961E-02 

-4.Q752E-02 

-7.19676-02 

-6.25146-02 

-5.49CU-02 

-3.12376-02 

-2.0456F-C2 

-2.5039E-O2 

-9.42366-02 

-8.20516-02 

- 7.50776-02 

-9.42366-02 

-8.20516-02 

- 7. 50776-02 

-5.52226-02 

-5. 75C0E-02 

-5.8C836-Q2 

-5.S919E-02 

-5. 01 306-02 

-5.86576-C2 


(b) Output. 


lit i 


Piqure i 0 . — Concluded. 

V.-ky ; 

-0 ; * • 


t TRAN IN W4ACH-1 . , MOPT-I , TAUB- . 1 ,TAUW».04 , XMTB- . 5 , XMTW- . 5 , ANGLE-4 5 . , 
SSMAX- . 321S , XRLE- . 25 , TR» , 4 , CRT- . 4 , XLBASE- .86 .XbOUTP- . 05 , 6 END 


( 4 ) Input. 


CAlCULAUCN UF SURFACE AND FLCw MfcLO PRESSURE D I S !R 1 8uT IONS 
F»7k FLCW A f FREE STREAM MACH NUMBERS Al CJH NEAR CNE . BELOW 
THE LO-fcR CRITICAL. LR ABOVE THE UPPkk CRITICAL ABCVJT A FINITE 
THICKNESS fcir.C-IUCSNTFC CIRCULAR ROCY COMBINATION WITH THE 
fcOC ! VAl t NT BODY OF RCVCUiTKh EITHER USER-SPECIFIED GR HAVING 
OKOlNAIcS * PROPORTIONAL TO X/L-(X/LI**N OR l-X/L-il- X/L !••* . 
THE RING HAVING A CONSTANT TH ICKNESS/CHCRC RAIIC, TAPER RATIO 
BcTWttN 0 AND l, AND WITH QNO I NATES l PROPORTIONAL TC X0AR/C 
-|XPAP/CI**P l« L-XeAR/C-ll-XBAR/C I **m IV USINC THE TRANSONIC 
EOMlVALtNCt RULE AND THE LOCAL LINEARIZATION *€THCC 


WING -BODY COMBINATION Cf.CHETRY ANC Fit* FtcLC CHARACTERISTICS 


EQUIVALENT RCLY THICKNESS ratio - 
EQUIVALENT BODY HAXIRU*- THICKNESS AT X/L - 
EXPONENT N FOR EQUIVALENT BCOY CROTNAIkS » 
S60* * I X I « 0 AT X/L « 

WING MAX. THICKNESS AT XBAR/C ■ 

WING THICKNESS/CHUMC RATIC » 

EXPONENT M FOR WING OkClNATFS ■ 

LEADING EDGE CF WING RCUT CHCRC AT X/l • 
TRAILING EOGE CF WING ROOT CHCRC AT X/L » 
PLANFCRW TAPER RATIO » 

LEA01NG EDGE PIERCES BCCY AT X/L » 

TRAILING EOGE PIERCES eUCY AT X/L - 
80QY RASE AT X/L • 

LEAOING EOGE SWEEP ANGLE IDfcGl - 
TRAILING EOGE SWEE« ANGLE I DEG I • 

LCCATICN OF WINGTIP LEADING ECCE AT X/L « 
LOCATION Of WINGTIP TRAILING EDGE AT X/L * 
NORMALIZED MAX. SEMISPAN SSMAX/L * 

ANGLE CF ATTACK ALPHA I0ECI =• 

RATIC CF SPECIFIC HEATS • 

FREE STREAM MACH NUMBER » 


C.1000CE 00 
C.50000E 00 
C.200G0E Oi 
0.21 1326 00 
C.5C000F. 00 
C • 40000b -0 1 
0.20000E 01 
C.2500CE 00 
C • 55000c 00 
C.40000E 00 
C.2S125E 00 
C. 571566 00 
C • 66QQQE 00 
C.4500QE 02 
C. 237556 C 2 
C.5715CE 00 
C.6915C6 00 
C.3215CE 00 
C.COOOO 
C.1400U6 01 
C.10000E 01 


START CF INTtGRATICN FRCM St 0 • ' ( X I • C TO NOSE 


X/L R BOCV/L THfc T A (DEG I CPIBCCYI CPIR/O- 1.00ICPIA/U> 2.00ICPTK/U- 3.001CPIR/C- 4.001CPIR/D- S.OOICPIR/D- 6.001 


0.2113 

0.2113 

0.03)3 

0.0333 

c.ucoo 

5.CCC0E 

01 

2.1 307E-C2 
2. 1307E-C2 

3.3155E-C2 

3.31596-02 

3.42 7CE-02 
3.42706-02 

3.4476E-02 

3.4476E-02 

3 .45486-02 
3. 45486-02 

3.45026-02 
3. 45826-02 

3.4600E-02 

3.46006-02 

0.2003 

0.2003 

0. C 320 
0. C32U 

O.GCOO 
9 . CCOOE 

01 

2.BBI jE- 02 
2. Edl 3E-02 

J.B1 71c-02 
3.81716-02 

3. 71246-02 
3. 7124E-02 

3.O065E-02 
3.6 0696-02 

3.52476-02 

3.52476-02 

).45e?E-C2 

J.4587E-C2 

3.40396-02 

3.40396-02 

0.1503 

0.1503 

0.0255 

0.0255 

c.ocoo 

9.0CO0E 

01 

6 • 669 2 E-C2 
6.6692E-02 

5.5471c -02 
5 • 94 7 IE -32 

4. 74 725- C2 
4. 74 726-02 

4.O071E-O2 

4.0071E-02 

3 .4755E-C 2 
3.47556-02 

3.C613E-C2 

3.U613E-02 

2 .72216-02 
2.T22U-02 

0.1C03 

0.1003 

C.0181 

0*0181 

O.OCOC 
<3 .0000c 

01 

1.13746-01 

1.1374E-C1 

7.56726-02 

7.56726-02 

5.0P74E-02 
5. OU 746-02 

3 .6 U0E-O2 
3 .6 i 1 BE -02 

2.56C7E-C2 

2.5oC7E-02 

1.74416-02 
1 . 744 IE-02 

1.07656- 02 

1 .07656- 02 

0.050) 

0.0503 

0. CO 56 
J .0096 

C.UCOO 

9.UOQOC 

01 

i.enec-ci 

1 . 6 1 34E-01 

7.94646-02 

7.94646-02 

4.01 33E-02 
4.01336-02 

1 . 70 38E -02 
1.7038E-02 

6.3612E-C4 
6 • 3612 E-04 

-1.205C6-C2 
-1 .2090E-02 

-2,24506-02 

-2.2490E-02 

0.004 3 
0 .OOt 3 

0.00C9 

0.0009 

o.ocoo 

9.100QE 

01 

3. 7747 E- Cl 
3. 7747E-01 

. 4.61076-02 
4.61076-02 

-7.91 176-03 
-7.91 176-0) 

-3.95126-02 

-3.5512E-02 

-6.19326-02 
-6 . 19 326 -02 

-7.53236-02 
- 7 • 93236-02 

-9.3533E-G2 

-9.3533E-02 

START CF 

INTTGRATICN VRCR 

S£A* 

MX) « 0 TC 

tail 






X/L ‘ 

KBOUV/L 

ThM A 1 DEG 1 

CPI bCCYl 

CPfR/O- 1.001CPIR/D- 2.0U1CP1R/0- 3.001 

CP1R/D* 4 .00 ICP ( R/D» 5.C0ICPIR/0- 6.001 

0.2113 

0.21M 

O.C33) 
0. 0 > )3 

C.uCO-J 

N.OCOUt 

01 

C. 13076-02 
2.13076-02 

J.3155L-02 

3.1159E-C2 

3.42706-02 
A.42 7CE-C2 

3.44766-02 

3.4476c-02 

3.4548E-0* 
3. 4548E-02 

3.45626-C2 
3 . 45 62E-C2 

J.4600C -02 
3.46006-02 

0.2503 

0.250) 

O.C375 
3.0 75 

C.OCOO 

5.00006 

01 

-2.4545E-C1 

-3.4545E-C3 

1.45506-02 
1 .49406-02 

2.30306-02 
2. 50306-02 

2.73116-02 

2.73116-02 

3.02786-02 

3.02766-02 

3. 23556-02 
3.2559E-C2 

3.44 13E-02 
3.44136-02 

0.3003 

0.300 

C.C420 

0.C420 

O.OCOO 

9.0COOE 

01 

6.461 7£-Cc 
-7.84j5t-02 

2.29116-02 

-1.196JE-02 

1.57346-02 
1.11 16E-02 

2.37666-02 
1 .99366-02 

2.74816-0 2 
2.5327E-02 

3.C616E-C2 

2.5237E-C2 

) . 32 78F-C2 
3.23206-02 

0.3503 

0.3503 

0.C440 
C. 0440 

O.OCOO 
9 . QCOOt 

01 

-;.4om-ci 
-t. 32186-01 

l. 00006 06 
-6.55726-02 

-4.44fcdE-0i 

-2.4B6AE-02 

l .9U5CE-Q3 
-6.79346-03 

5 . 4154k- 03 
4. 53586-03 

1 . 5926c- C2 
1. 28036-02 

2 . 13CSC-C2 
L • 9335E-02 

0.4003 

0.4003 

0.04 34 
0.04 34 

C.OCOO 

9.00006 

01 

- 1 . 655 36-C1 
-1.41646-01 

-1 • 2i2b6-Cl 
-9. 04926-02 

- £. 02 62E-C3 
-4.67576- 02 

-9.66516-03 
-2 .46106-02 

-2.26251-03 

-1.05046-02 

< 5.01C36-C3 
-2.38666-04 

1.14706-02 

1.03186-03 


Figure 11.- Sample input/output for a wing-body combination having a 
circular body and with TR f 0, 6 t# > 0, o ■ 0. 


Ill 


0.4SU1 

C.C413 

C.OCQC 


- 1,335 K-Ql 

-1.93506-01 

l. 00 001 04 

-1.31% 16-02 


-3.41786-01 

2.60 IK-01 

0.4901 

0.041} 

9.00006 

01 

-S.961M-02 

-I. 24322*0/ 

-9.20/76-02 

-1.23446-02 

-1,03 106-02 

-9.54161*03 

-i.iicat-oi 

0.5003 

0.C391 

o.ocoo 


- 7.4*556- (2 

•1.43 I2€-Ci 

-1./94U-01 

-V .6QT06-Q2 

-l.'IPlbC-Oa 

-1. 46481*0/ 

-5.*!C7t-ei 

0*5003 

0.C391 

S.O0OOI 

oi 

-l,54.»6?-02 

-3.19/56-0/ 

-4.68511-02 

-3.55656-0/ 

-/.593/f-O? 

-1.01406-0/ 

-l.l?19f-C/ 


Sf Aft r Cf 4U*C*IS0MC CUCULAUC* 

&UKRICNK tntuvMKfc turn at x/t • 0.54/3/6 do 

*A &4O0Y/1 THffAlUfcC) CPfeOtV) CPI9/0- l.OOKHB/O- 2,OOUPIt»/n* 3*tl>K>»ll»/i> 4.QCKMH/D* 5.00 fCPIf»/C» 6.00I 


0.5501 

0.1301 

O.CS90 

0.0390 

0.0000 
4 .00006 

Oi 

9.15696-0) 
7. WUt-02 

-1./M6F-01. 
-1, 56546-02 

-l. 47516-01 
-4.54/7002 

1 .00006 06 
-4.4273C-02 

- ■*• 994 76- C2 
-1.85276-02 

- i* 45 806-02 
-1.21C66-02 

-3.11446-02 

-2.64921-02 

C ,6CC j 
0.4003 

C . C4 1 0 

0.0414 

o . tc co 

4.00006 

Oi 

*.t5l5f-C2 

1.14446-02 

1 .94140-01 
1.57466-02 

-1.38631-01 

-1*26926-02 

-2 .17146-01 
-2.65916-02 

-U06C26-01 
-3*1 3196*02 

- 7 . t6 4 56-CJ 
-3.16676-02 

-6 • 105*6-02 
-1.03936-0? 

0.4503 

0.4901 

0.04)9 

0.0439 

O.OOJO 

9 .ooooe 

01 

r/.«77*6-Cl 
-2. 77*71-01 

-i .uote-oi 

,•1.7*271-01 

3.94P8C-02 
-1 • 15526-01 

-4.2 i 966-01 
-a .*69 fit-02 

-o . VI 996-02 
-7.60596-02 

-0.O8716-C2 
-6 . 70046-02 

-5.47796-0/ 

-0,04406-0/ 

0*7003 

0.700) 

0.0420 

0.0420 

o.ocoo 

«,0000£ 

01 

-1 .31051-01 
-l .31096-01 

-1 .07766-01 
-1.0T76C-01 

-4.25376-02 
-9,253 76-02 

-8.39716-02 

-8.35716-02 

- 7.79506 -C2 
-7.79506-C2 

- 1 » 12476-02 
-I.I257F-Q2 

“4.99C2F-C2 

-6.45C/O02 

0.7501 

0.7903 

0.0315 

0.03J5 

C.UCOJ 
9, OCOOf 

01 

-1.16056-01 
-1 .16096-01 

-9.772600/ 

-5.77246-0/ 

-6.5Jslt-C/ 
-4, 975| 6-02 

-8.55726-02 

-8.54726-02 

-0.26496-C2 

-8.26456-02 

-0.0406-02. 

-8.0406-02 

-*«.0376t-O; 

-7.43780-0/ 

0,4CU3 

0.4003 

4.'. €320 
O.C3/0 

o.ocoo 

9.0C00E 

01 

-0.01.46-02 

-8.«l/4t-C2 

-I.99//6-0/ 

-7.39221-3/ 

-e.oc <,76-04; 
-4.00576-0/ 

-8.1 /« 700/ 
-8.12276-02 

-S.2K26-02 

-8.2IC/6-07 

-4.28C«f-G2 

-6.28C46-C2 

-4* 37866. -02 
- 8. 13866-02 

0.4903 

0.0903 

0,0255 

0,0255 

o.ocoo 

9 . OCCQt 

01 

- 3.95356-C2 
- i.9» J5C-02 

-4 • 704C6-02 
70406-02 

-5.919100/. 

-5.901002 

-t> ,t>t»8C6 -02 
-6 .06806-02 

- 7 . 205 76“ C2 
-7.20576-02 

-7.6247E-C2 

-7.62476-02 

?96rC< 
- 1.96796-0/ 


C«AC CCtMICIEKT • C.IC 44 C' OC 


(b) Output. 


Figure U.- Concluded. 

ns - 



faTRANIN AMACH= 1 . , MOPT= 1 , TAUB= . I , TAUW« . 04 , XKTB** . 5 ,XMTU° . 5 ,MiGLE«4S . . 
SSMAX* .336 , XRLE- . 25.TR= .4 .CRT- .4 , X LEASE" . 86 , XLOUTP- . 05 , NL-3 . , fa END 


(a) input. 


CALCULATION CF SURFACE AND FLOW FIELD PPtSSL«E Ci SIR! BUT IONS 
FUR Fli:» AT FREE STREAM m AC H NUMBERS AT l)R Kt AH ONE, BELCH 
THE LChEK CvniCAl, OH A Ho V E THE UPPER CRITICAL AB CU T A FINITE 
THlCKAtSS WlNO-INCENTEC mucv curuinatipn with the bocy having 
an elliptic cross Section that maintains a ccnstaki ratio of 

MAjCfl/»INCR AXES ALONG THE ENTIRE BCOV LLNJTH WITH THE EQUI- 
VALENT MOOT OF REVOLUTION t ITrtc* USER-SPECIFIED OR HAVING GH- 
01 NATES' R PRUPQRl loNAL TO A /L - ( X/ L 1 • *N UR l - X/L - I i - X/L I • *N , 

THE m|NG HAVING A CONST AM THICKNE SS/CMCRC HATIC, TAPER RATIO 
Hi T n E El* 0 AND i, ANO WITH JRCINATES l PROPORTIONAL TO XtJAR/f 
- ( X PAR /C I •• m OR l-xeAR/C-U-XPAR/CT**M f>Y USING THE TRANSONIC 
EQUIVALENCE RULE ARC THE LOCAL LINEARIZATION *£TMCC 


WlNG-BCOY COMB IN AT ION OECMcTHV AND FLCW FltLC CHARACTERISTICS 


RATIO OF SEMImAJOR/SEMIMINCR AXIS * 
EQUIVALENT BODY THICKNESS RATIO « 
EQUIVALENT BOOT MAXIMUM THICKNESS AT X/L * 
EXPONENT N FJR EQUIVALENT OCOY C'pUINATES * 

see* 'ix> • o at x/l 

M 1 NG MAX. THICKNESS AT X0AR/C » 

WING ThICKNESS/CHCRO RATIO * 

EXPONENT M FOR WING CRflNATES * 

LEADING EDGE OF RING RCOT CHCRl AT X/l = 
TRAILING EDGE CF WING MOOT CHORC AT X/L » 
PlANFCRR TAPER RATIO « 

LEAOINC EDGE PIECES 8CCY AT X/L - 
TRAILING EUGE PIERCES ECCY AT X/L * 

BODY BASE AT X/L » 

LEAOING fcDGt SWEEP ANGLE ICcG) =• 

TRAILING E06t SWEEP ANGLE IOEGI « 

LOCATION OF WINGMP LEAOING EDGE AT X/L « 
LOCATION OF WINGTIP TRAILING ECCE AT X/L « 
NORMALIZED MAX, Sf: M I SPAN SSMAX/L - 
ANGLE CF ATTACK ALPHA ICEGI * 

RATIO OF SPECIFIC HEATS - 
FREE STREAM MACH NUMBER “ 


C.3000UE 01 
C.IOOOCE 00 
C. 500001 00 
C.20000E 01 
C .21 1326 00 
C.50C0CE 00 
C • 400GCE - 0 1 
C.200C0E 01 
C.250CCE 00 
C.S5OC0E 00 
C.4000CE 00 
C.22612E 00 
C. 58919c 00 
C.96QCCE 00 
C.450C0E 0/ 
C.25054E C2 
C.58800E 00 
C.7080CE 00 
C.338CCE CC 
C.COOOO 
C.14000E 01 
C.10000E 01 


START CF 

INTEGRA 

TICN FPCM 

see* 

•m * c tc 

NCSE 






X/L 

RBOCY/L 

THE TAtPEGI 

CPItlCCY » 

CPIR/C- l.OOICPTR/O* 2 • 00 ) CP |H /D s 3.0U1CPIR/C* 4.00ICPIR/C* 5.00ICPIH/C- 6.001 

0.2113 

0.2113 

0.0577 

0.0192 

C.JCOO 

S.OCCUt 

01 

4.79 7tr.-C2 
1. 2s* 6E-C2 

3. 7658E-02 
J.025JE-C2 

3.5267F-C2 
3. 3 360t-02 

3.4V21E-02 

3.sO5)E-02 

3.47976-02 

3.43C6E-02 

3.474CE-02 

3.44266-02 

3.47 10E- 02 
3.44 9 IE-02 

0.2003 

0*2003 

C. C55 5 
0.0185 

C.CCCC 
9 .OOOJf 

01 

5 . 71CSE-C2 
l » 27« 7 E- 02 

4 .2837E-02 
1.50206-02 

J.81U8E-02 

3.6l60t-02 

3.6535E-02 

3.50236-02 

3.55C0E-C2 
3. 49S36-02 

3.47«4£-C2 

3.4423E-02 

3.41546-02 

1,39296-02 

0.1503 

0.1503 

0.C442 

O.Cls7 

o.ococ 

s.OOOOE 

ul 

1 .C31 )t-Cl 
5 . 09t. 1 E- 02 

6. J955E-02 
5 . 5 d 796-02 

S.8544L-02 
4.6461 E-02 

4 »054 3t-02 
3. 9oilE-02 

3 . 50 206-02 
1.4495E-02 

3.07026-02 

3.0446E-02 

2.7338F-02- 
2. 71046-02 

0.1003 

0.1003 

0.0313 

o.oic*. 

J.UGOQ 

9.UCO0E 

01 

1 ,595dE-01 
1. .s,:06-C2 

7.67936-02 
7 .2 BdOc -02 

5. to3sfc-02 
5. 01346-02 

3.O454C-02 

3.5786E-02 

2.5796E-02 

2.54196-02 

1.75626-02 

1.7321E-C2 

1.00486-02 
1 . 0681 E-02 

0.050 3 
0.0503 

0.C1 c6 
U. CC*>5 

o.occo 

9.0C00E 

Oi 

1. 2 7B6E- Cl 
1 . 51oOE-Ol 

e. 06116-02 

7.8353C-02 

4.C410E-02 
3 . 9851 E- 02 

l .7164E-01 
l .6912E-02 

7.07116-04 
5. 65246-04 

-1.2045E-02 

-1.2136E-02 

-2.2459E-C2 

-2.25226-02 

0 . 00 s J 

0.00N3 

0. C01 r > 
U. C0C5 

O.OCOC 

9.0C00E 

01 

4.445J6-C1 

3.4004E-01 

4,6ll5r.-C2 
4 .60966 -02 

- 7. SCddE-03 

- 7. 9147E-03 

-3.961CE-C2 

-3.9513E-02 

-6.1932E-02 
-6 • 19 33E-02 

-7.53236-02 

-7.53246-02 

-9.35326-C2 

-9.3533E-02 

START CF 

/NTi'GRA 

I If N FMCM 

SE8‘ 

MX> * 0 1C 

1 AIL 






A/L 

R'H OV/L 

THt TAIOtGI 

CMIBCCYI 

CP(w/D» I.OUICPIH/C* 2 . 001 CP Ik/U* 3.C0ILPIR/D- 4.0CICPIR/0- 5.C0ICP(H/C« 6.001 

0.2113 

0*2113 

C. C‘a t 7 

o.ois? 

C.CCCC 

S.OCCOc 

01 

s . 79 »4£-0*. 
1.24L8C-C2 

3 . 1 6 5oe - 02 
3.02536-02 

3. 52676-02 
3. 136UE-Q2 

3.s92iE-02 
3 .405 36 -02 

3. 4 7976-02 
3.4 J06E-02 

3.474CE-C2 

3.S426E-Q2 

3.471 Of -02 
3 .44 S IE - 02 

0,2503 

0.2503 

0. L«5C 
0. 0/ 1 T 

c .ocoo 

5.0CC0E 

01 

1. 75-.3E-C2 
-8. 6543E-0) 

1 • 6 3266 -02 
1.32&S6-G2 

2.3756E-C2 

2.24136-02 

2.7627E-02 

2.7017C-02 

3.04 55E-C2 
J.0109E-02 

3.26716-C2 

3.2449F-02 

3.44S2E-C2 

3.4337E-02 

0.30C3 
O.i 003 

C. 0728 

0.024) 

C.OLOO 
9.CC0 )r 

ul 

- 1 . 58 19F-02 

- 3. 29* CE-OZ 

- 7.59356-03 
-B. 08926-03 

7 . l287t-03 
6. 9l2sk-OJ 

1 .5o62E-02 
1 »p 55 26 -02 

2.1651E-C2 
2 . 1 62 7E-02 

2.L359E-C2 

2.63176-02 

3.0170E-C2 
3 .01406-02 

0.3503 

0.3503 

0.0784 

C.C261 

o.tcoc 

9.0C0JE 

01 

-t. I022f>02 
-l.l 39« f-Ul 

1.C0006 06 
-4.73256-02 

1. 30796-03 
-1.49106-02 

7.18486-03 
i • 2 75 )E -04 

1.36041-02 

9.6654E-C3 

1.91406-02 

1.66166-02 

2.3AWE-C2 

2.2133E-02 


Figure IJ.- Sample Input/output (or a wing-body combination having an 
elliptical body and *ith ?R + 0, B ta > 0, <» - 0. 
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0.-.003 

0*400 ) 

J.079J 

0.0? 

t .0100 
>*.00001 

01 

-i.e)02t-ct 
-l .6? >»f-C* 

-1.4C64? -01 
- 0 . * l ».i -02 

- 1. 704.-02 
-5. 

-l .6 tC9t-C2 
-1.0H»«I5-G.» 

- 7.20546-03 
-i. *. 8551-0? 

1.30274-2* 
- »*54 » «f -13 

8.7*.e7F-C.* 

5* 4*40, -0 1 

0*<.50 ) 
0*4*503 

o. 1 ft: 

j.n?54 

O.JCOO 

s.ocooc 

Oi 

- 4 . 

-l . N 29f-Ul 

- 1.6 7 4 71 -0. 
‘*45l -72 

i • Old > Of. 
• '..247/:- J2 

-4 .•.)4,n«jr-o? 

-l.tSU! -0. 
-2.4/ 7*1-02 

- 7. *24. %-ii * 

-1.401 ip-c? 

- * . . s V * - J~ 

- *.o‘-**l* -0 1 

0*5003 

0*5C03 

j.o7u 

J.02 37 

C.OCLC 
9. 00006 

01 

- f . 30 fct :-c; 

-4.45CM - 0 ? 

-1. : C e .* 1 - J1 
-5.619** -02 

-l.?4. ;:-c* 
- * . /O49r-02 

-4.124.' :-j.* 

- J 734t -Oi 

-1.142. t-Ci 

-2.*»945» - 02 

-1 .24 IAI-C7 
-i. 73156-0? 

- J.pfc 12* -« ■ 
-*.u;7 - 0 ? 


S t *Q f CF StPfaUJMC CALCLL4Ili\ 

513*6 0 5 'Ml C*lCU«4fUN Sl4«li 41 w/t * 0 ,isl}2t JO 


*/i oinur/L r hi i a tutu i w^i -««. cv » cpir/d* 1 .ook pi 1 */. • * .joicpif /o* j.oimMrf/c* *.ou i c p c o/r.« « .00 k pi*/: • *.<*•>» 


0.550 J 
0.550 J 

J. 07.62 
O.C227 

0«>l COO 

o.occ.o* 

01 

- .. 7c Ml-Oj 
s.*»4?4F-Ci- 

-4. ?d«* 7-. - J? 
-2.1445* -0. 

-1.761 *F-Ol 
- *. .. 7if -C*. 

1 . JOUOE 06 
— 3 . HOC IF —02 

- 1.49 3/1 "0. 

- 31 43t-C/ 

- J. 1% J56-C2 
-2.6C7T*.-C2 

- 2 .'. Wr.-C2 

16r -C. 

0.60C3 
0.600 1 

C.07C7 
3. C« 36 

Q.’ICUO 
5. ’KOOt 

0* 

l.c6>*<il-l| 
4.21 »J* -0. 

i. 71141-0. 
5 . /-. 74| 

— / .40 hOC— 01 
5.*. 7 7 4 f-U « 

-i.»C64fc-01 
-l .45776-02 

-1 .*01 U-Cl 

-..4471L-C2 

-■*. 343 7f-C. 
-2.71 336- C 2 

;■ -*;< 
-2.7110- - 2 

O.ftSCl 

0.6503 

0.0738 
C.C2 46 

c .ocoo 

4 • C COOt 

J» 

-1.4-051 -01 
- 1 .ftSftbf-Cl 

- i . J/0. s :. - J 2 
-1 .04**1! -Jl 

.. n-,o*-02 
-**.0l->2f-02 

-2.1 4446-01 
-6 .4 336E-02 

- 1 . 404 5t — *)/ 
-ft . 2H9.*f -02 

-7.40645-0/ 
-5. 74955-07 

-ft «**27Ct -C2 
-5. 3fe4»* -02 

0.7003 
0.7 0 j 3 

O.C72* 
J. c?4i 

O.'JCOQ 

4.0COI)t 

01 

-4. 7124: -U. 
— 4 . 4 1 <*».* -«:i 

-3. *.7/4* - 01 
- 3. :•»/«• -Cl 

-2. *>**436-01 

9. *-60/6-02 

— 1 .‘*4 lit -01 

1.U35CE-C/ 
- 1 .2 34 TF- C 4 

-2. 3764E-C2 
-1.05246-01 

- 3.6C<*1 -C? 
-S.2«;9* -T2 

O.?«03 

0.7403 

J. Co47 
7.0216 

c.ocoo 
4 . 0 C OOt 

01 

-5.46C7F-1. 
-i.Jl 34t * 01 

- • .4 - ■} } 

•I.W ') f - C 2 
-4.041 4 F- J2 

-6.525 it-C 2 
-d .5 56 ^F -0? 

-**.24716-02 

-6.?62'Jt-02 

-6.02646-02 

-P.05136-02 

- 7. 

- 7 • H# 54* -u. 

0.6003 

0.4003 

0.C554 

0.01 a? 

Q.OCOw 
4.0 0006 

04 

-5.47361-02 
-4. *21 ?r-C2 

- 7 .4/461 -02 
-t.?0rt*.c-02 

-7.50 Jlt-02 

- t. 1064*- 02 

-6 .0 7oC6-U« 
-3.1674* -02 

-6.ie4it-ci 

-6.2357C-C2 

-6.26176-C2 

-6.2968C-02 

-6.327C -C2 
- 92 

0.650 J 
0.850? 

7.0441 

J.oitf 

C.OCOO 

4.30001 

01 

-2. c 8 l>* -0 7 
-5.7 

-•*.2523* -02 
- 6.062 *1. -02 

-7.6121 F-02 
-0. 02 CO f-02 

-ft .».?0 r »t -02 

-6 .7 1 34c -02 

-f.*753t-02 
-7.2117c -02 

-7.60766-02 

-7.64146-02 

-7.45t2» -C2 
-7.4744*-02 


CRAG CCtHICIfM • 0. *«S1TE-U1 


(b) Output. 


**»i«ui** Concl'ideJ. 
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.... space activities of' the United States shall be,, 
ondtteted so as to contribute to the expansion of human knowt* 
edge of phenomena in the atmosphere and space. The Aimimitraiiad j 
shall provide for the widest practicable and appropriate dissemination, 
pf odor matt on concerning its activities and the results thereof 

—National Aeronautics and Space Act of 1958 
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lilll 


TECHNICAL NOTES: Information less broad 
in scope hm nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
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